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ABSTRACT 

Initial  yield  surfaces  of  unidirectional  fiber  reinforced 
composites  subjected  to  combined  longitudinal,  transverse 
normal,  and  longitudinal  shear  loads  are  calculated.   The 
composite  is  composed  of  filaments  arranged  in  a  periodic 
square  array  embedded  in  a  matrix  material.   The  constituent 
materials  are  assumed  to  be  homogeneous  and  isotropic,  while 
the  composite  is  assumed  to  be  macroscopically  homogeneous 
and  transversely  isotropic.   The  finite  element  method, 
using  linear  strain  triangles,  is  employed  to  calculate  the 
stresses  throughout  the  composite.   Von  Mises  yield  criterion 
is  used  to  calculate  the  elastic  limit  of  the  local  micro- 
scopic combined  stresses  in  the  composite.   A  parametric 
evaluation  is  carried  out  by  changing  individual  constituent 
properties  and  evaluating  the  effect  of  this  variation  on 
the  composite  properties.   Yield  surfaces  for  several 
functional  composites  are  included,  as  well  as  a  computer 
program,  in  FORTRAN  IV  language,  for  the  calculation  of  a 
yield  surface. 
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I.   INTRODUCTION 

A.   COMPOSITE  MATERIALS 

The  requirements  of  modern  technology,  especially 
that  of  the  aerospace  industry,  have  necessitated  the 
development  of  materials  possessing  much  higher  strength 
to  weight  ratios  than  common  structural  materials.   This 
result  has  been  achieved  through  the  introduction  of  com- 
posite materials,  that  is  two  or  more  constituent  materials 
combined  in  such  a  way  that  the  gross  properties  of  the 
composite  are  favorable  combinations  of  those  of  the 
constituents.   This  indicates  that  the  desirable  properties 
such  as  strength/weight  and  stiff ness/weight  ratios  are 
increased  while  less  desirable  properties  such  as  .low 
ductility  are  minimized. 

A  specific  type  of  a  two  constituent  composite  is  a 
unidirectional  fiber  reinforced  composite.   This  consists 
of  continuous  lengths  of  a  stiff  reinforcing  fiber  imbedded 
in  a  ductile  matrix  material.   All  the  reinforcing  fibers 
are  oriented  in  the  same  direction  within  a  single  layer  of 
material.   Generally  several  layers,  or  laminae,  are  bonded 
one  on  top  of  another  to  form  a  high  strength  laminated 
composite. 

Each  of  the  lamina  may  be  considered  as  an  anisotropic, 
homogeneous  material.   The  transverse  properties,  i.e., 
those  in  a  direction  perpendicular  to  the  filament  orienta- 
tion, are  weaker  than  those  in  the  longitudinal  (fiber) 
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direction  for  a  single  lamina.   To  counteract  this  weakness 
in  a  laminated  composite,  adjacent  lamina  are  arranged  such 
that  they  have  different  filament  orientation.   The  funda- 
mental unit  is  therefore  the  single  layer  of  unidirectionally 
reinforced  materials.   It  is  this  type  of  basic  layer  that 
is  analyzed  in  this  investigation. 

The  physical  dimensions  of  such  a  lamina  can  vary 
considerably,  although  the  largest  is  generally  but  a  fraction 
of  an  inch  thick.   The  fibers  vary  in  diameter  from  .0001" 
for  glass  and  .0003"  for  graphite  to   .004"  for  boron.   Lamina 
thickness  depends  on  the  materials  used  and  the  intended 
purpose.   An  individual  lamina  can  be  produced  in  nearly  any 
desired  thickness.   However  standard  laminae  are  becoming 
popular.   For  example,  most  graphite  comes  in  a  standard  .01" 
thickness.   Allowing  for  the  presence  of  the  matrix  material 
such  a  lamina  may  contain  in  excess  of  25  fibers  across  its 
thickness. 

B.   DEFINITIONS 

The  terminology  of  composite  materials  is  often  confusing 
due  to  a  lack  of  precision  in  defining  specific  terms.   For 
example, at  times  it  is  unclear  whether  a  term  such  as  tensile 
strength  refers  to  the  ultimate  strength  of  a  material,  i.e., 
its  failure  point,  or  to  the  point  of  initial  yield.   In 
composite  materials  it  is  necessary  to  be  especially  clear 
by  noting  whether  a  term  refers  to  a  composite  property  or 
to  a  property  of  a  constituent  material.   In  order  to  avoid 
such  difficulties  the  definitions  of  this  section  will  be 

adhered  to  in  this  study. 
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1.  Macroscopic  Quantities 

These  are  quantities  of  the  composite  material. 
They  are  calculated  by  assuming  that  the  composite  consists 
of  anisotropic,  homogeneous  layers. 

2.  Microscopic  Quantities 

These  quantities  are  determined  on  a  point  by 

point  basis  throughout  the  composite.   The  different 

mechanical  properties  of  the  matrix  and  filament  materials 

are  considered  in  the  calculation  of  these  quantities.   To 

see  the  relationship  between  macroscopic  and  microscopic 

quantities  consider  Fig.  (1).   Figure  (la)  represents  a 

body  with  microstress,  a  ,.   In  Fig.  (lb)  the  macrostress 

S   in  the  x-direction  obtained  from  the  microstress  field 
x 

of  Fig.  (la)  is 

a 

S   =  -  |  °x(y)  dy  (I_1) 

o 

S   is  the  macrostress  for  the  face  x  =  b,  and  represents 

the  average  value  of  the  microstresses . 

3.  Initial  Yield  Strength 

That  combination  of  macroscopic  stresses  at  which 
a  point  in  the  composite  first  reaches  its  elastic  limit, 
i.e.,  no  plastic  strain  has  occurred  in  the  composite. 

This  study  is  restricted  to  plane  macrostress 
fields  with  components  Sv,  S  ,  and  Sv„,  where  the  z  axis 
is  along  the  length  of  the  fiber.   Specific  cases  of 
initial  yield  strength  associated  with  this  investigation 
are 
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a    (b,y) 


>     X 


a)      Microstress   Field   a    (b,y) 

A 


O  b 

b)   Macrostress  S 
Figure  (1) .   Relationship  of  Microstresses  and  Macrostresses 

a)  Longitudinal  Yield  Strength,  when  the  only 

macroscopic  stress  is  in  the  axial  direction,  i.e., 

S  ?    0,  S   =  S   =  0. 
Z      '   X     xz 

b)  Transverse  Yield  Strength,  when  the  only  macro- 
scopic stress  is  in  the  transverse  normal  direction,  i.e., 

S   7*  0,  S   =  S    =0. 
x  '   '   z    xz 

c)  Longitudinal  Shear  Yield  Strength,  when  the 

only  macroscopic  stress  is  in  the  longitudinal  shear 

direction,  i.e.,S   7*0,  S   =S   =0. 

xz       x    z 
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4 .   Second  Deviated  Stress  Invariant ,  J 

The  second  deviated  stress  invariant  at  a  point 
is  defined  as   [lj 


-   I  (axy2  +  V2  +  °zx2)  +|  [(",   "  VS 


(oy  -  az)2  +  (a,  -  o„)2|  (1-2) 


This  quantity  can  be  related  to  the  yield  strength  of  a 
material . 

5.   Effective  Micros tres s ,  o 

The  effective  microstress ,  o,  is  proportional  to 
the  square  root  of  the  second  deviated  stress  invariant,  [l] 


a   -  I   [(o-a 
/2  L   x 


y)2  +  (a  -  az)2+  (oz  -  ox)2 


+  6  (oxy2  +  oyz2  +  ozx2) 


(1-3) 


6.  Specific  Strength 

The  ratio  of  a  particular  initial  strength  term, 
as  defined  above,  to  the  specific  density  of  the  composite. 
The  strength  considered  must  be  specified,  for  example, 
specific  longitudinal  strength.   This  may  be  applied  to 
composite  materials  as  well  as  single  phase  materials. 

7 .  Specific  Modulus 

The  ratio  of  the  Young's  modulus  of  a  material 
to  the  specific  density  of  the  material.   In  the  case  of 
anisotropic  materials  it  must  be  specified  which  modulus 
is  being  considered,  such  as  the  specific  transverse  modulus 
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8.  Ultimate  Strength 

That  point  in  macroscopic  stress  space  at  which 
either  a  matrix  element  or  a  filament  element  actually 
fails,  i.e.,  breaks, cracks ,  or  separates.   Due  to  the 
brittle  nature  of  filament  materials  the  initial  yield 
strength  and  ultimate  strength  usually  coincide  for  filament 
materials. 

9.  Initial  Yield  Surface 

In  a  general  three-dimensional  system  the  stress 

state  is  defined  by  six  stress  components,  three  normal 

stresses,  S  ,  S  ,  and  S  .and  three  shearing  stresses,  Sv„, 
x    y        z  3  xy ' 

S   ,  and  S   .   Thus  an  initial  yield  surface  may  be  repre- 
sented in  macrostress  space.   Any  macroscopic  stress  state 
within  the  surface  means  that  the  effective  microstress  at 
every  point  in  the  composite  is  less  than  that  required  for 
yield.   Any  stress  state  lying  on  the  surface  means  that 
some  point  in  the  material  has  an  effective  microstress 
equal  to  the  elastic  limit.   For  brittle  materials  without 
a  plastic  region  loading  states  outside  this  surface  are  not 
possible.   For  ductile  materials  stress  states  outside  the 
initial  yield  surface  (within  the  limits  of  the  ultimate 
strength  of  the  material)  will  result  in  plastic  deformation 

The  three-dimensional  yield  surface  considered  in 
this  study  is  a  specific  case  of  the  general  initial  yield 
surface  described  above.   Macroscopic  stresses,  S  ,  S   ,  and 

S..„  are  set  equal  to  zero  in  the  calculation  of  the  yield 
y  » 

surface.   The  resultant  surface  represents  the  yield  surface 
due  to  plane  macrostress  loading  conditions. 
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10.   Stress  Concentration  Factor 

A  stress  concentration  factor,  denoted  as  SCF, 
is  generally  defined  as  the  maximum  microstress  in  the 
composite  to  the  applied  macrostress.   In  this  study  three 
particular  cases  of  stress  concentration  factors  are 
investigated. 

a)   SCFX.   Transverse  loading  stress  concentration 
factor,  associated  with  macrostress  state  (Sx,  S  ,  Sxz)  = 


(Sx,  0,  0). 


(CJx)max 

SCFX  =  2^£  (1-4) 

S 
x 


b)   SCFZ.   Longitudinal  loading  stress  concentration 
factor,  associated  with  macrostress  state  (Sx,  S  ,  Sxz)  = 


(0,  S  .  0)  . 

v   '  max 
SCFZ  = (I_5) 

z 

■ 

c)   SCFXZ .   Longitudinal  shear  loading  stress 
concentration  factor,  associated  with  macrostress  state 


<Sx'  Sz'  Sxz>  =  <0'  °<  Sxz>« 


(°xz) 

s 


SCFXZ  =  ^^—  (1-6) 


xz 

11.   Transverse  Stiffness,  ET 

The  transverse  stiffness  of  a  composite  is  the 
macroscopic  modulus  of  elasticity  in  a  direction  perpen- 
dicular to  the  fiber  direction. 
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12.  Longitudinal  Stiffness ,  E 

L 

The  longitudinal  stiffness  of  a  composite  is  the 
macroscopic  modulus  of  elasticity  in  the  fiber  direction. 

13 .  Longitudinal  Poisson's  Ratio ,  vT 

L 

The  longitudinal  Poisson's  ratio  is  the  measure  of 


the  contraction  in  the  transverse  direction  due  to  tension 

EX 

e 


in  the  longitudinal  direction,  i.e.,  v   =  —   for  t^e  case 

(sx,  sy,  sz)  =  (0,  o,  sz). 

14.  Transverse  Poisson's  Ratio,  v 

1 T 

The  transverse  Poisson's  ratio  is  the  measure  of 

the  contraction  in  the  transverse  direction  due  to  tension 

in  the  same  plane,  that  is  v   =  e  /e  ,  for  the  case 

l     y   x 

(Sx'  Sy'  Sz)=  (Sx'  °'  0)  • 

15.  Shear  Modulus,  G 

The  composite  shear  modulus  is  the  modulus  of 
elasticity  in  shear  due  to  longitudinal  shear  loading.   As 
the  transverse  composite  shear  modulus  is  a  dependent 
quantity  any  reference  to  composite  shear  modulus  in  this 
study  will  mean  the  longitudinal  shear  modulus  unless 
specifically  stated  to  the  contrary. 

16 .  Filament  Volume  Fraction,  v^ 

The  filament  volume  fraction  is  the  ratio  of 

filament  material  volume  to  total  material  volume, 

v   -   Filament  Volume 
f      Total  Volume  (1-7) 

Note  that   0  <_  vf  <_  1.0. 
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17.   Matrix  Volume  Fraction,  v 
m 

The  matrix  volume  fraction  is  the  ratio  of  the 
matrix  material  volume  to  the  total  material  volume, 


Matrix  Volume 
Total  Volume 


m       Tnfal  Vnlnmp  (1-8) 


Again,  0  <_   v  <_   1.0.   For  a  material  without  voids  the  sum 
of  the  matrix  and  filament  volume  fractions  must  equal 
unity, 

vf  +  vm  =  1.0.  (1-9) 


C.   BACKGROUND 

When  a  tensile  test  is  performed  on  a  composite  material 
specimen  to  determine  an  ultimate  tensile  strength  the 
value  obtained  is  a  macroscopic  property.   When  design 
specifications  require  certain  properties  these  also  are 
macroscopic  properties.   To  accurately  determine  the  behavior 
of  a  composite  and  to  bridge  the  gap  between  design  require- 
ments and  design  practice,  macroscopic  behavior  must  be 
related  to  microscopic  behavior  within  the  composite. 

Numerous  analytical  approaches  for  the  determination  of 
composite  properties  have  been  utilized.   These  include 
strength-of -materials,  stress  function,  finite  difference, 
and  finite  element  techniques,  among  others.   Chamis  and 
Sendeckyj  [2]  present  a  comprehensive  review  of  the  various 
analytical  methods. 

Specific  examples  of  micro-stress  analysis  include  the 
longitudinal  shear  loading  solutions  of  Adams  and  Doner,  [3] 
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and  Tsai ,  Adams  and  Doner   [4].  The  longitudinal  load 
case  was  considered  by  Bloom  and  Wilson  [5]  using  a  series 
solution  technique.   Transverse  normal  loading  was  investi- 
gated by  Tsai,  Adams  and  Doner  [4]  ,  and  Adams  and  Doner  [6] 
using  finite  difference  analyses,  and  by  Foye  [7]  using  a 
finite  element  method.   Lin,  Salinas  and  Ito  [8]  employed 
finite  element  techniques  to  consider  the  problem  of 
combined  loading,  by  superposing  the  effects  of  longitudinal, 
transverse,  and  longitudinal  shear  loading.   Lin  et  al 
calculated  an  initial  yield  surface  based  on  this  combined 
loading  analysis  and  later  extended  their  work  into  the 
plastic  region,  [9],[lo].   Adams   [Lij  also  considered  inelastic 
loading  but  restricted  his  study  to  transverse  normal  loading 

The  combined  loading  problem  has  received  the  least 
attention  though  it  is  the  most  practical  from  a  design  stand- 
point.  It  is  of  great  importance  to  the  designer  to  know 
what  combination  of  macrostresses  will  result  in  a  micro- 
stress  distribution  sufficient  to  cause  yield.   The  initial 
yield  surface  provides  this  information. 

D.   SCOPE  AND  OBJECTIVES 

The  analytical  approach  used  in  this  investigation  was 
originally  developed  and  programmed  by  Lin,  Salinas,  and 
Ito.  [8]   The  author's  objectives  have  been 

i)   Modify  the  program  to  allow  use  of  an  overlay 

system  structure  in  order  to  reduce  the  required  computer 
storage. 

ii)   Substitute  the  use  of  random  access  disk  unit  for 

the  original  tape  storage  methods. 
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iii)   Perform  a  parametric  study  of  the  effects  of 
constituent  properties  on  the  properties  of  the  composite. 
iv)   Compare  several  realistic  material  combinations  to 
evaluate  their  structural  efficiency. 

v)   Provide  selected  initial  yield  surface  plots  for 
design  assistance. 

vi)   Make  some  general  observations  relating  composite 
behavior  to  constituent  materials. 
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II.   METHOD  OF  ANALYSIS 

The  analytical  method  described  in  this  section  was 
developed  by  Lin,  Salinas,  and  I to  [8].  The  present  discus- 
sion is  a  general  presentation  of  the  basic  principles  of 
that  analysis.   Appendix  A  is  a  detailed  description  of 
the  method. 

The  basic  assumptions  necessary  to  the  analysis  are, 
i)   The  filament  and  matrix  materials  are  isotropic 

and  homogeneous. 

ii)   The  filament  and  matrix  materials  are  perfectly 

bonded.  • 

It  should  also  be  noted  that  the  macroscopic  stress 

state  is  one  of  plane  stress,  with  components  S  .  S„,  and 

x    ^ 

Sxz* 

The  microscopic  stress  distribution  within  a  unidirec- 

tionally  reinforced  composite  depends  upon  the  arrangement 
of  the  filament  material  within  the  matrix.   In  order  to 
facilitate  this  analysis  the  filament  array  is  assumed  to 
be  a  doubly  periodic  rectangular  array.   In  such  an  array 
the  two-dimensional  (x-y)  surface  may  be  divided  into 
identical  segments  called  basic  blocks,  of  width  2a  and 
height  2b  as  shown  in  Fig.  (2).   Away  from  the  edges  of  the 
composite  the  microscopic  stress  distribution  should  be 
identical  within  each  block.   Due  to  symmetry  within  a  basic 
block  only  one  quadrant  need  be  analyzed,  as  shown  in  Fig.  (3) 
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Figure  (2).   Rectangular  Array  of  Filaments  and  a  Basic  Block  [8 J 

y 


Figure  (3).   A  Quadrant  of  the  Basic  Block  [o] 
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In  order  that  no  voids  be  formed  and  that  adjacent 

blocks  neither  separate  or  overlap  it  is  necessary  that 

the  rectangular  basic  blocks  remain  rectangular  during 

deformation.   In  order  to  produce  a  combined  macroscopic 

plane  stress  state,  as  well  as  maintaining  the  rectangular 

form  of  the  blocks,  four  specific  load  cases  are  defined. 

These  four  cases,  shown  in  Fig.  (4),  are  superposed  to 

give  plane  stress  and  boundary  compatibility. 

Load  cases  I  and  II  are  plane  strain  states  with 

e   =0.   Load  case  III  is  a  generalized  plane  strain  case 

with  unit  strain.   Load  cases  I,  II,  and  III  are  superposed 

to  give  S   =0.   Load  case  IV  is  a  state  of  longitudinal 

shear.   All  four  load  cases  have  Sxy  =  Syz  =  0.   The  boundary 

conditions  for  each  problem  are  shown  in  Fig.  (4). 

These  four  load  cases  are  solved  by  the  finite  element 

method  using  linear  strain  triangles.   Felippa  [l2]  presents 

a  detailed  discussion  of  finite  element  theory  and  a 

derivation  of  the  linear  strain  triangle.   The  finite 

element  model  utilized  for  the  circular  filament  cross 

section  problems  is  shown  in  Fig.  (5).   This  model  was 

adapted  from  that  of  Lin  et  al   [9].  In  this  work  convergence 

tests  were  conducted  to  establish  the  validity  of  the  model. 

Each  load  case  yields  a  microstress  distribution  and 

associated  macrostresses .   The  microstresses  are  designated 

by  0°^  •  ,  with  subscripts  determining  the  stress  component 
13  I 

and  superscripts  denoting  the  load  case.   For  example  o ^ 

is  the  microstress  a±j  at  a  specific  point  associated  with 
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Figure  (5) . 


Finice  Element  Model  for  Circular  Filament 
Cross  Section. 
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problem  I.   Macroscopic  stresses  are  represented  by  S^ 
with  subscripts  and  superscripts  similar  to  microstresses . 

Appendix  A  treats  the  solution  of  the  individual 
problems  and  demonstrates  how  the  results  may  be  combined 
to  relate  the  combined  loading  microstresses  at  a  point  to 
macrostresses  as 

[o]  =  [A]  [S]  ("-1) 

See  Appendix  A,  Eq.  (A-13) ,  for  the  explicit  form  of  these 

matrices. 

AS  noted,  the  microstresses  vary  from  point  to  point. 

Therefore  the  A  matrix,  evaluated  as 

La]  =  M  [»]"'  (II"2) 

is  also  point  dependent. 

in  this  thesis  the  Von  Mises  yield  criteria  was  taken 
to  govern  initial  yield  of  the  composite.   In  terms  of 
microscopic  stresses,  yielding  occurs  when  [l] 


j  =  a   2  +  a   2+  a   2+  - 
2   xz     ¥      y    3 


V+V  +  tf 


a  J  +a„z  +  a/  -  ax 


ay  "aYaz  "azaxJ 


1  Y  2  (H-3) 

3   © 


where  Y   is  the  elastic  limit  of  the  material  for  the  tensile 


o 
stress  state 


Using  Eq.  CH-2)  to  relate  microscopic  and  macroscoprc 
stresses  Appendix  A  shows  that  the  yield  equation,  Eq.  (II-3) , 


may  be  written  as 
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ClSz2  +  C2SxSz  +  C3Sx  +  C4Sxz2  =  I  Yo2  (II-4) 

The  terms  C  ,  C~,    Co,  and  C.  are  second  order  polynomials 
of  the  components  of  the  A  matrix  and  are  functions  of 
position.   For  given  values  of  S   and  S   ,  Eq.  (II-4)  may 
be  solved  for  the  value  of  S   necessary  to  cause  yield. 
Each  point  within  the  fundamental  block  will  yield  a  dif- 
ferent value  of  S  .   Each  combination  of  S   and  Sv_  yields 
two  values  of  S   at  each  point,  a  maximum  and  a  minimum. 
The  minimum  of  the  maximum  values  and  the  maximum  of  the 
minimum  values  determine  two  yield  points  for  the  given 

values  of  S„  and  S   .   Continuing  in  this  manner  the  initial 
x      xz  3 

yield  surface  for  a  composite  subjected  to  a  macroscopic 
plane  stress  state  is  calculated. 
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III.   EVALUATION  OF  NUMERICAL  RESULTS 

A.   INITIAL  YIELD  SURFACE  CHARACTERISTICS 

The  two-dimensional  plots  of  the  yield  loci  of  this 

investigation  have  the  general  form  of  Fig.  (6) .   Each 

locus  represents  a  specific  value  of  longitudinal  shear 

stress,  Sv„,  with  the  outermost  locus  corresponding  to 

Sv_  =  0,  as  seen  in  Fig.  (6). 

Each  point  on  a  yield  locus  indicates  that  some  point 

within  the  material  would  yield  under  the  given  macrostress 

state.   The  analysis  used  in  this  investigation  provides 

the  location  of  the  area  of  yield  associated  with  each 

point  of  a  yield  locus.   These  area  locations  refer  to  the 

triangular  elements  of  the  finite  element  grid,  Fig.  (5). 

Figure  (6)  is  the  yield  surface  for  a  50%  Boron/6061 

Aluminum  composite.   The  elemental  areas  that  first  reach 

yield  under  specific  loading  conditions  are  listed  in  Table 

I.   The  corresponding  macrostress  states  are  noted  on  Fig. 

(6) .   Proceeding  in  this  manner  a  particular  triangular 

element  could  be  associated  with  every  point  on  each  yield 

locus.   This  information  could  be  useful  to  a  designer  in 

determining  if  failure  is  associated  with  a  particular 

region. 

On  initial  consideration  it  seems  somewhat  unusual 

that  a  yield  locus  is  not  symmetric  about  the  S   ancj  g   axes 

£       z 

To  explain  this  behavior  first  recall  the  yield  equation  of 
the  previous  section,  Eq.  (II-4) . 
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C!Sz2+   C2SzSx   +   C3Sx2    +   C4Sxz2    =   \  V 


(III-l) 


The  exact  definition  of  the  terms  C-,  ,  C2/  C~,  and  C.  can  be 
ascertained  from  Eq.  (A-22) . 


TABLE  I.   INITIAL  YIELD  FOR  S    =  0 

xz 


Point 

S  (ksi) 

2% 

S  (ksi) 
z 

Yield  Location 
(Element) 

A 

0 

114 

66 

A' 

0 

-114 

66 

B 

+  15 

+  112 

67 

B' 

+  15 

+  75 

88 

C 

+  30 

+  100 

67 

C 

+  30 

+  .45 

88 

As  stated,  the  points  of  the  initial  yield  surface  are 
calculated  by  setting  Sxz  and  Sx  to  specified  values  and 
calculating  the  value  of  Sz  resulting  in  a  microscopic  stress 
sufficient  to  cause  yield.   Assume  that  the  macrostress 
combination  (Sxz/  Sx,  Sz-,)  is  known  to  be  such  a  yield  point. 
It  will  be  shown  that  for  different  combinations  of  S   and 

A 

S  ,,  i.e.,  (Sx,  -  S  -i )  etc.,  a  specific  type  of  symmetry 


zl 
exists . 


For  (Sx,  Szl)  Eq.  (III-l)  becomes 


C.S^t2  +  C0S„iS„  +  CoS  *  +  CiS   *=  i  Y  2 


'l°zl 


'2°zl°x 


'3°x 


■4°xz 


2_  1 


(IH-2) 
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For  (-Sx,  -Szl)  Eq.  (III-D  is 

ClSzl2  +  C2SzlSx  +  C3Sx2  +  C4Sxz2=  J  Yo2       (IH-3) 

Equations  (III-2)  and  (III-3)  are  identical  and  show 
that  yield  occurs  at  the  same  point  and  for  the  same 
absolute  values  of  macrostresses  in  quadrants  I  and  III 
of  the  yield  locus,  as  shown  in  Fig.  (7).   This  indicates 
symmetry  about  the  S    axis,  perpendicular  to  the  (S  ,S7) 
plane. 

For  combinations  of  (Sv,  S  ,)  such  that  their  product  is 

X      Z  J. 

negative,  (Sx,  -S  ^)  or  (-S  ,  S  ,  ) ,  Eq.  (III-l)  becomes 

clSzl2  -  C2SzlSx  +  C3Sx2  +  C4Sxz2  =  i  YQ2       (III-4) 

Comparing  Eqs.  (III-2)  and  (III-4)  it  is  seen  that  if  (Sxz, 
Sx,  Szi)  is  a  valid  solution  of  Eq.  (III-l)  as  was  assumed, 
then  (Sxz,  S  ,  ~SZ^)  and  (S  z,  -S  ,  Sz^)  cannot  be  a  solution. 
For  a  positive  S   and  negative  Sz,  solution  of  Eq.  (III-l) 
will  yield  a  different  minimum  S  ,  say  Sz2-   This  corresponds 
to  a  different  point  in  the  basic  block  than  does  the 
solution  for  quadrants  I  and  III.   The  combinations  (Sx,  -Sz2/ 
Sxz)  and  (~SX,  Sz2,  Sxz)  identically  satisfy  Eq.  (III-l) 
establishing  the  equality  of  yield  points  in  quadrants  II 
and  IV  of  Fig.  (7).   The  symmetry  of  the  yield  locus  about 
the  Sxz  axis  is  thus  established. 

Another  interesting  phenomenon  apparent  in  some  yield 
surfaces  is  the  occurrence  of  square  ends  on  the  yield  loci 
corresponding  to  small  values  of  Sxz.   This  is  shown  in  Fig. 
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(8) .   Generally  this  indicates  that  yield  is  occurring  in 
the  filament  material  rather  than  the  matrix.   Note  that 
the  S   scale  is  considerably  larger  than  the  Sv  scale. 
The  longitudinal  stress  values,  Sz,  in  the  vicinity  of  the 
square  ends  are  many  times  larger  than  the  transverse,  Sv, 
values.   The  Sz  values  are  so  dominant  that  increasing  the 
transverse  stress  from  point  A  to  point  B  of  Fig.  (8)  has 
little  effect  on  the  yield  surface.   Once  the  macrostress 
state  B  is  reached  yield  occurs  in  the  matrix  material. 


Figure  (7) .   Symmetric  Nature  of  a  Yield  Locus 

As  stated  the  yield  locus  plot  is  a  two-dimensional 

representation  of  a  three-dimensional  surface.   Thus  the 

separation  between  successive  yield  loci  is  an  indication 

of  how  steep  the  slope  of  the  yield  surface  is  in  that  area. 

The  dark  areas  of  the  plots  indicate  an  extremely  steep  slope, 

such  that  adjacent  loci  actually  plot  very  close  to  one 

another. 
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B.   EVALUATION  OF  COMPOSITE  MECHANICAL  PROPERTIES 

1.  Composite  Strength 

Composite  strength  properties  were  determined 
directly  from  the  yield  surface  associated  with  each 
particular  composite. 

i)   Longitudinal  yield  strength  was  taken  as 
points  S   on  the  initial  yield  surface  when  S   =  Sx=  0. 

ii)   Transverse  yield  strength  was  taken  as  points 
on  the  initial  yield  surface  when  Sz  =  SxZ  =  0. 

iii)   Longitudinal  shear  yield  strength  was  taken  as 
the  point  S   ,  on  the  initial  yield  surface  when  Sx=  Sz  =  0. 

2.  Composite  Elastic  Constants 

Composite  elastic  constants  were  calculated  from  the 
results  of  the  individual  longitudinal,  transverse,  and 
longitudinal  shear  loading  problems.   The  macroscopic 
stresses  associated  with  each  loading  case  were  used  in 
evaluating  the  transversely  isotropic  stress-strain  relations 
developed  in  Appendix  B.   As  discussed  in  this  appendix  a 
transversely  isotropic  material  is  an  orthotropic  material 
with  a  single  plane  of  isotropy.   It  should  be  noted  that 
the  assumption  of  transverse  isotropy  limits  the  applicability 
of  this  method  of  determining  the  elastic  constants  of  a 
composite.   Composites  are  not  generally  transversely 
isotropic.   The  particular  case  treated  in  this  investigation 
is  that  of  a  square  array  of  filaments.   For  this  particular 
fiber  distribution  the  composite  properties  approximate  those 
of  a  transversely  isotropic  material.   For  more  general 
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filament  arrangements,  such  as  a  rectangular  array,  further 
consideration  must  be  given  to  the  applicability  of  these 
stress-strain  relations.   The  analytical  method  employed 
to  calculate  the  displacements,  stresses  and  strains  is 
not  effected  by  this  restriction  however.   The  calculation 
of  these  quantities  is  applicable  to  a  general  rectangular 
filament  array. 

For  convenience  the  transversely  isotropic  stress-strain 
relations  are  repeated  here. 


a)   e* =  i;  (sx  -  Vy>  -  fJ  s2 

b)     £y  =  jL  (sy  -  Vx)   -    3J  Sz 

C>       £z    =    -   F1    <Sx    +    Sy)    +   J-  Sz  ("1-5) 

L  ^L 

XY        G      bxy 

T         x 


e)       £yz   =  _L_  s 


GL      Y* 


-  J-  S 

'L 


xz         GT       xz 


In  the  transverse  loading  case,  problem  I,  the  macro- 
scopic strain,  ex,  corresponding  to  the  unit  displacement 
imposed,  was  calculated.  The  macroscopic  strains  e  and 
e   are  equal  to  zero. 

In  the  longitudinal  loading  case,  problem  III,  the  macro- 
scopic strain  cz   is  equal  to  unity.   In  this  problem 

Ex  =  £y  =  0. 
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Equations  (III-5) ,  evaluated  for  these  two  loading 

conditions,  were  solved  for  longitudinal  and  transverse 

composite  stiffness  and  Poisson's  ratios,  E  ,  E  ,  v  ,  v  . 

Li    T    L    1 

The  longitudinal  Poisson's  ratio  was  evaluated  directly 

from  Eq.  (III-5c)  in  terms  of  load  case  I  as 

I 


v. 


S1*  S1 
x   y 


(III-6) 


Longitudinal  stiffness  was  evaluated  from  Eq.  (III-5c)  in 
terms  of  load  case  III  as 


111      /„IH     IIIv 
E   =  S7   -  vT  (S„    +  St      ) 

L     z 


(IH-7) 


Lx"x     wy 
Equation  (III-5b)  may  be  evaluated  in  terms  of  load  cases 


I  and  III  and  the  resulting  equations  solved  simultaneously 
give 


E, 


E   = 
T 


V   =  VT 
T     L 


(III-8) 


(III-9) 


The  longitudinal  composite  shear  modulus,  G,  was  obtained 
from  the  results  of  the  longitudinal  shear  loading  case, 
problem  IV.   The  shear  strain  corresponding  to  the  imposed 
unit  displacement  was  calculated  and  used  in  conjunction 
with  the  longitudinal  shear  strain  Eq.  (III-5f ) . 
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3 .  Specific  Strength  and  Specific  Stiffness 
Specific  strength  and  specific  stiffness  were 

calculated  by  dividing  the  respective  properties  of  the 
composite  by  the  composite  density.   Densities  of  the 
various  material  combinations  were  calculated  from  the 
specific  weight  data  of  Table  II,  which  lists  properties 
for  all  materials  used  in  this  study. 

4 .  Stress  Concentration  Factor 

Three  specific  values  of  the  stress  concentration 
factor  are  calculated  in  this  investigation,  corresponding 
to  three  specific  loading  cases.   All  comply  with  the 
general  definition  of  the  ratio  of  the  maximum  microscopic 
stress  to  the  macroscopic  stress. 

i)  SCFX  refers  to  the  transverse  loading  problem, 

problem  I.   It  is  defined  as  the  maximum  microscopic  stress, 

av,  divided  by  the  average  surface  traction,  S  .   From  Eq. 

(1-4), 

(ax) 

SCFX  =     max  (111-10) 

x 

ii)   SCFZ  refers  to  the  longitudinal  loading  problem, 
problem  III.   It  is  defined  as  the  maximum  microscopic 
stress,  a  ,  divided  by  the  average  surface  traction.   From 

Eq.  (1-5), 

(CTz) 

crip?  max 

SCFZ (III-ll) 

bZ 
iii)   SCFXZ  refers  to  the  longitudinal  shear  loading 
problem,  problem  IV.   It  is  defined  as  the  maximum  micro- 
scopic shear  stress  av_,  divided  by  S„_,  the  average  shear 
stress  of  this  problem. 

39 


(axz> 

max 
SCFXZ  = (111-12) 

xz 


C.   COMPARISON  OF  NUMERICAL  RESULTS 

In  order  to  establish  the  validity  of  the  analysis 
used  in  this  investigation  comparisons  were  made  with  avail- 
able theoretical  and  experimental  results.   It  is  not  the 
intention  here  to  establish  the  superiority  of  one  method 
of  analysis  over  another  but  rather  to  ensure  that  the 
solutions  are  compatible. 

1.   Transverse  Stiffness ,  E 

Using  a  finite  difference  approach  Tsai  et  al  [4  J 
determined  the  transverse  stiffness  of  a  40%  fiber  volume 
Boron/6061  Aluminum  composite  as  Em  =  20  x  106  psi.   The 
present  analysis  gives  E   =  20.97  x  106  psi. 

Adams  and  Doner,  [p]    also  using  a  finite  difference 
analysis,  calculated  a  theoretical  value  of  transverse 
stiffness  for  a  47%  E-glass/Epoxy  model  as  1.55  x  106  psi, 
and  reported  an  experimental  value  of  1.30  x  106  psi.   This 
analysis  gives  1.42  x  10e  psi.   Adams  and  Doner  also  list 
a  comparison  between  their  results  for  Glass/Epoxy  and 
several  other  analyses,  including  finite  element,  complex 
variables,  and  other  solution  techniques.   These  results 
are  compared  to  the  present  analyses  in  Table  III. 

Transverse  stiffness  for  Boron/Epoxy  and  Glass/Epoxy 
also  fell  within  the  variation  bounds  calculated  by  Dow 
et  al  J15]  .   Further  comparisons  were  made  with  results  of 
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TABLE  II.   MATERIAL  PROPERTIES 


Specific  Wt 
Material1  E(psi) v     Y^(psi)2   (lb/in3) 


1. 

Boron 

60xl06 

0.2 

4.8xl05 

.085 

2. 

S-glass 

12xl06 

0.2 

5.0xl05 

.090 

3. 

E-glass 

10.6xl06 

0.22 

4.0x10s 

.090 

4. 

MODMOR  I 
(graphite) 

60xl06 

0.25 

2.5xl05 

.072 

5. 

MODMOR  II 
(graphite) 

40xl06 

0.25 

4.0xl05 

.063 

6. 

THORNEL  25 
(graphite) 

25xl06 

0.25 

1.8xl05 

.052 

7. 

THORNEL  40 
(graphite) 

40xl06 

0.25 

2.5xl05 

.052 

8. 

6061  Aluminum 

lOxlO6 

0.3 

.346xl05 

.100 

9. 

2024  Aluminum 

8.1xl06 

0.34 

.130x10s 

.100 

10. 

NARMCO  23  87 
(epoxy) 

.5xl06 

0.3 

.300x10s 

.045 

11. 

4617  Epoxy 

.78xl06 

0.34 

.800xl05 

.045 

1.  Materials  1  through  7  are  filament  materials. 
Materials  8  through  11  are  matrix  materials. 

2.  Y   is  the  elastic  limit  of  the  material  for  the 

o 

tensile  stress  state. 
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Ekvall  |l6]  ,  Whitney  |i.7]  and  others.   The  majority  of 
these  results  are  in  the  form  of  plotted  curves  which 
makes  direct  numerical  comparisons  difficult  but  does 
provide  a  general  comparison  of  results. 

TABLE  III.   COMPARISON  OF  TRANSVERSE  NORMAL  LOADING  SOLUTIONS 


Ef/Em=20 
Analysis 


vf='2    vm='34 


vf=.55 


Transverse  Stiffness,  Er 


Adams  and  Doner  [_6j 

Wilson  (l3] 

Pickett  (l4] 

Foye  [7] 

Present  analysis  (47%) 

Present  analysis  (60%) 

Present  analysis  (50%  estimate) 


1.86xl06  psi 

1.86xl06  psi 

1.85xl06  psi 

1.83xl06  psi 

1.42xl06  psi 

2.05xlOG  psi 

1.81xl06  psi 


2.   Composite  Shear  Modulus 

Adams  and  Doner  [3] ,  using  a  finite  difference  method, 
have  calculated  numerous  theoretical  values  of  composite 
shear  moduli  for  specific  material  combinations.   For 
Glass/Epoxy  they  obtained  a  shear  modulus  of  .417xlOG  psi 
and  .589xl06  psi  for  40%  and  55%  models  respectively.   The 
present  analysis  yields  .447xl06  psi  for  a  47%  model  of  the 
same  material. 
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For  Boron/Epoxy  Adams  and  Doner  obtained  .  441xl06  psi 
and  .658xl06  psi  for  40%  and  50%  models,  compared  to  the 
present  values  of  .476xl06  psi  for  a  47%  model. 

Direct  numerical  comparisons  were  available  from 
Blakslee  et  al  [l8]  who  utilized  Whitney's  micromechanics 
model  [l7]  based  on  an  Airy  stress  function.   The  particular 
graphite  investigated  was  "Thornel"  (a  commerical  graphite 
produced  by  Union  Carbide) .   For  Thornel  25  in  a  50%  model 
Blakslee  reported  a  theoretical  value  of  the  composite  shear 
modulus  of  .518xl06  psi  and  an  experimental  value  of  .60xl06 
psi.   This  compares  to  a  present  value  of  .531xlOe  psi. 
Similar  values  for  a  60%  model  of  Thornel  40  are  .931xl06  psi 
theoretical  and  . 740x10 6  psi  experimental,  compared  to  a 
present  value  of  .768xl06  psi.   These  comparative  values  for 
Thornel  are  summarized  in  Table  IV. 

Graphical  comparisons  were  available  with  Ekvall  [16]  , 
Tsai  et  al,  [4]  and  Whitney   [L7]  . 

TABLE  IV.   SHEAR  MODULUS  COMPARISON  FOR  THORNEL/EPOXY 

Shear  Modulus,  G(psi  x  10~6) 
Analysis THORNEL  25*     THORNEL  40** 


Blakslee  [l8] 

.518 

.931 

Experimental  [l8] 

.600 

.740 

Present  Analysis 

.531 

.768 

*  THORNEL  25   Ef/Em=  32    vf=.25    Vm=.34    vf=.5 

THORNEL  40   E^/E  =  51    vf=.25    V  =.34    V4r=.65 
t   m         t         m         t 


** 
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The  graphical  results  of  Fig.  (9)  were  compiled 
by  Reidinger  et  al  [l9]  .   This  compares  the  values  derived 
by  Tsai  [20]  ,  Dow  et  al  [21]  ,  and  Ekvall  [16]  for  Glass/Epoxy 
The  results  of  this  analysis  are  marked  by  X.   Both  Tsai  and 
Dow  et  al  obtained  their  results  from  consideration  of 
variational  bounds.   Tsai  however  employs  a  "contiguity" 
factor  dependent  on  the  amount  of  contact  between  filaments. 
The  present  analysis  and  that  of  Dow  et  al  assume  that  the 
filament  is  completely  surrounded  by  matrix  material.   This 
accounts  for  the  closer  correspondence  of  the  present  results 
to  that  of  Dow  et  al. 


2.0 


•H 

a 

I 

o 

X 


o 


.H 
O 

£ 
u 
o 

CO 


1.5 


1.0 


Tsai  [2  0] 

Dow  et  al  [l5]         J1 

Ekvall  [16]  /  / 

x   Present  Analysis  // 


0.5 


0   .1   .2   .3   .4   .5   .6   .7   .8 
Fiber  Volume  Fraction  (v^) 


E  =10.6xlOb  psi   Ern=0.5xlCb  psi  vf=.22  vm=.35 
Figure  (9).   Shear  Modulus,  G,  vs.  Fiber  Volume  Fraction,vf 
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3.   Longitudinal  Stiffness ,  ET 

L 

Longitudinal  loading  has  not  been  considered  as 
frequently  as  transverse  or  shear  loading,  resulting  in 
fewer  numerical  comparisons.   Blakslee  J18]  again  gives 
theoretical  and  experimental  values  for  Thornel/Epoxy 
composites.   For  the  50%  volume  Thornel  25,  he  reports  a 
theoretical  value  of  12.78  x  10 e  psi,  and  an  experimental 
value  of  11.07  x  10 6  psi.   This  compares  to  a  value  of 
11.96  x  106  psi  in  this  study.   For  65%  Thornel  40  similar 
results  are  theoretical  26.20  x  106  psi,  experimental 
24.6  x  106  psi  and  a  present  value  of  24.35  x  106  psi.   This 
data  is  summarized  in  Table  V. 

Reasonable  correlation  was  also  obtained  in 
comparison  with  graphical  results  of  Tsai  [2  0j  and  Ekvall  [l6] 

TABLE  V.   LONGITUDINAL  STIFFNESS  COMPARISON  FOR  THORNEL/EPOXY 

Longitudinal  Stiffness,  EL(psi  x  10~6) 
Analysis THORNEL  25* THORNEL  4  0** 

Blakslee  [l8]  12.78  26.20 

Experimental  [l8]  11.07  24.60 

Present  Analysis  12.44  24.35 

*  THORNEL  25   E./E  =32    v. =.25    V  =.34    v^=.5 

f '  m       r        m        r 

**  THORNEL  40   Ef/Em=51    Vf=.25    Vm=.34    vf=.65 
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4 .   Composite  Yield  Strength 

Strength  comparisons  were  more  difficult  to  obtain. 
In  some  cases  strength  quantities  described  by  identical 
terminology  and  applied  to  seemingly  similar  materials  dif- 
fered by  several  hundred  per  cent.   One  reason  for  such 
discrepancies  is  the  lack  of  consistent  definition  of 
strength  terms.   Valid  comparisons  only  have  meaning  if 
the  quantities  investigated  are  precisely  defined. 

Another  important  point  is  the  effect  of  the  matrix 
yield  strength  on  composite  yield  strength.   Matrix  mate- 
rials of  the  same  class,  such  as  different  brands  of  epoxy, 
normally  have  elastic  constants  that  show  little  variation 
from  one  another.   The  difference  in  yield  strength  however 
can  be  considerable.   The  result  of  this  difference  is  shown 
in  Table  VI,  comparing  a  40%  filament  model  of  Boron/2024 
Aluminum  with  a  matrix  yield  of  1.30  x  10**  psi  and  Boron/6061 
Aluminum  with  a  matrix  yield  of  3.46  x  101*  psi. 


TABLE  VI.   EFFECT  OF  MATRIX  YIELD  STRENGTH 
ON  COMPOSITE  PROPERTIES 

v  =.  5 

Boron/6061     Boron/2024 
Aluminum Aluminum %  Difference 

EL  (psi)    20.53x106      18.86x106  8. 

ET  (psi)    31.07x106      28.71x106  8. 

G   (psi)     6.33x106       5.28x106  19. 

S   (psi)    94.53xl03      38.9xl03  143. 

L 

ST  (psi)    28.50xl03      11.20x103  154. 
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As  seen  in  Table  VI  the  change  in  elastic  properties 
is  insignificant  in  comparison  to  the  change  in  strength. 

The  change  that  does  occur  in  the  elastic  properties  is  due 

to  minor  differences  in  the  elastic  properties  of  the  two 

matrix  materials.   The  reason  for  the  very  large  increase  in 

strength  may  be  seen  from  Fig.  (10) . 

0 

A 


J2   „ 
B 


o 


Figure  (10).   Strength  Comparison  of  Different  Alloy 
Matrix  Materials. 

The  slope  of  curve  A  represents  the  composite 

longitudinal  stiffness.  E  ,  assumed  to  be  identical  for 

L 

two  composites,  each  having  the  same  filament  material 
but  with  different  alloys  of  the  same  matrix  material. 
As  seen  in  Table  VI  the  values  of  ET  in  such  a  case  are 
quite  similar.   Curve  B  represents  the  stress-strain  curve 
of  one  matrix  material,  curve  C  that  of  the  other  matrix 
material.   They  are  assumed  to  have  approximately  the 
same  elastic  modulus  (as,  for  instance,  two  different 
aluminum  alloys)  but  with  different  yield  stress.   For 
the  case  of  matrix  B  the  yield  stress  will  be  reached  at 
a  composite  strain  e  .   The  resultant  composite  stress  is 

Ll 
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If  the  matrix  material  is  that  corresponding  to  curve 
C,  matrix  yield  will  occur  at  composite  strain  z    ,    resulting 

in  composite  stress  S   .   Thus  a  variation  in  matrix 

L2 
strength  alone  may  markedly  affect  the  composite  strength. 

No  comparison  of  strength  terms  should  be  attempted 
unless  the  complete  properties  of  both  constituents  are 
specified. 

The  example  of  Boron/2024  Aluminum  of  Table  VI  is 
directly  comparable  with  the  results  of  Adams  [ll].   For 
identical  material  properties  using  a  finite  element 
analysis  with  constant  strain  triangles  Adams  reported  a 
theoretical  transverse  yield  strength  of  11,200  psi, 
identical  with  the  results  of  this  analysis. 

Longitudinal  strength  may.  be  readily  analyzed  by 
strength-of -materials  methods.   Holister  and  Thomas  [22J 
developed  a  relationship  for  determining  the  ultimate 
tensile  strength  of  a  composite  undergoing  axial  tension 
by  considering  the  relation 

ST  =  afvf  +  amvm  (111-13) 

J_i 

where  v_  and  v   are  the  filament  and  matrix  volume  fractions, 
f      m 

respectively. 

For  the  case  of  axial  load  for  continuous  fibers,  as  in 
this  study,  the  strain  will  be  constant.   Holister  and 
Thomas  considered  the  particular  case  where  the  strain  is 
such  that  the  ultimate  stress  of  the  filament  is  attained 
while  the  matrix  material  is  still  in  the  elastic  range. 
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This  is  equivalent  to  assuming  a  stress-strain  relation 
corresponding  to  Fig.  (11).   Curve  A  represents  the 
composite  stress-strain  relation  and  curve  B  that  of  the 
matrix. 


ult 


m 


ult 


Figure  (11) . 


eult 

Stress-Strain  Diagram  for  Failure  Occurring 
in  the  Filament. 


If  the  matrix  remains  in  the  elastic  region,  the 

strain  z    , . .  for  ultimate  strength,  S  , , ,  must  be  less  than 
ult'  3     ult 

that  necessary  to  cause  initial  yield  in  the  matrix.   Under 
these  conditions 

eult  =  Yo/Ef  ■  °A  (111-14) 

Equations  (111-13)  and  (111-14)  along  with  the  relation 


Vr  +  Vm  =   1 

f    m 


(111-15) 


give 


Sult  =  eult  vfEf d-Em/Ef)  +  £ultEm     (111-16) 


Reidinger  et  al  reported  an  equivalent  form  of  this 
equation  as 


ult    f  of    m  ■ 


(111-17) 
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where  a  is  defined  as  the  stress  in  the  matrix  at  the 

ult 
ultimate  strain  of  the  filament.   This  is  shown  in  Fig.  (11). 

The  yield  strength  of  the  filament  is  YQ  . 

f 

It  must  be  repeated  that  these  relations  are  for  a 
composite  that  fails  while  the  matrix  is  still  elastic.   It 
should  be  pointed  out  however  that  Holister  pursues  the 
matter  further  to  consider  the  case  of  the  matrix  being  in 
the  plastic  range. 

For  the  more  common  case  of  the  matrix  material 
yielding  prior  to  the  filament  reaching  yield,  relations 
similar  to  (111-13)  and  (111-14)  may  be  utilized  to  find 
the  initial  yield  strength.   In  the  case  of  Fig.  (12)  the 
strain  sufficient  to  cause  initial  yield  for  longitudinal 
load,  S  ,  is  marked  by  e_ . 

L  -Li 

a 


m 


Figure  (12) 
For  this  case 


Stress-Strain  Diagram  for  Failure  Occurring 
in  the  Matrix. 


e  = 
L 


m 


m 


m 


m 


*_£_ 
E* 


(111-18) 


where  YQ   is  the  axial  yield  of  the  matrix.   Using  Eqs. 

m 
(111-13)  and  (111-18)  yields 


SL  =  vf°f*  +  vmYo 


(111-19) 


m 
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The  stress  in  the  filament  at  the  yield  strain  of  the 
matrix  is  a   *,    where 


Ef 

°f*    =  Yo    E~  (111-20) 

m   m 

The  results  for  longitudinal  strength  calculated  in 
this  study  differed  by  approximately  2%  from  those  found 
by  Eq.  (111-20) . 

5.   Composite  Poisson's  Ratio 

a.  Longitudinal  Poisson's  Ratio.   The  longitudinal 
Poisson's  ratio,  vT ,  as  calculated  from  the  transversely 
isotropic  stress-strain  relations,  Eq.  (III-5) ,  exhibits  a 
linear  behavior  with  respect  to  filament  volume  fraction. 
The  results  of  Fig.  (13)  for  E-glass/Epoxy  correspond  to 
those  of  a  law  of  mixtures  relation. 

vL  =  vfvf  +  vmvm  (111-21) 

This  is  in  accord  with  the  theoretical  results  of  Rosen  et 
al,   [23]  Whitney  and  Riley,   [17]  and  Tsai  [go]. 

Specific  experimental  results  were  available  from 
Blakslee  et  al  Q.8]  for  THORNEL/Epoxy .   This  data  also 
showed  adherence  to  a  rule  of  mixtures. 

b.  Transverse  Poisson's  Ratio.   The  comparison  of 
the  transverse  Poisson's  ratio,  vT,  with  other  theoretical 
studies  raised  several  interesting  points.   Values  found  in 
the  literature  ranged  from  0.7   [29]   to  considerably  less 
than  0.1.  J18]   Thus  the  literature  indicates  a  range  of  the 
composite  transverse  Poisson's  ratio  that  far  exceeds  that 
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Figure  (13).   Longitudinal  Poisson's  Ratio  vs.  Filament  Volume 
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of  the  Poisson's  ratios  of  the  constituents.   The  higher 

values  exceed  the  theoretical  maximum  of  0.5  associated 

with  isotropic  materials.   [24J 

In  this  investigation  using  the  transversely  isotropic 

stress-strain  relations,  Eqs.  (III-5) ,  results  were  obtained 

for  v  .   With  the  addition  of  a  small  amount  of  filament 
T 

material  the  composite  transverse  Poisson's  ratio  takes  on 
a  value  greater  than  the  Poisson's  ratio  of  either  constit- 
uent.  For  the  materials  considered  the  increase  for  a  30% 
filament  volume  was  approximately  22%  above  that  of  the 
matrix.   In  all  cases  the  value  of  v   remained  below  0.5. 

The  transverse  Poisson's  ratio  decreased  with  increasing 
filament  volume.   At  high  filament  volume  v   values  were 
less  than  the  Poisson's  ratio  of  either  constituent.   The 
maximum  decrease  noted  in  this  study  was  47%  below  the 
filament  Poisson's  ratio. 

To  investigate  the  cause  of  the  transverse  composite 
Poisson's  ratio  having  values  outside  the  range  of  its 
constituents  several  additional  computer  analyses  were 
made.   The  constituent  and  composite  properties  of  these 
analyses  are  listed  in  Table  VII. 

If  the  transverse  Poisson's  ratio  was  dependent  only 
on  the  Poisson's  ratios  of  the  constituents  and  the  filament 
volume,  the  results  of  analyses  1,  2,  and  3  below  would  be 
identical.   As  noted  in  the  previous  section,  the  longi- 
tudinal Poisson's  ratio,  v  ,  is  dependent  only  on  these 

Jj 

quantities.   This  fact  is  borne  out  by  the  data  of  Table  VII. 
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TABLE  VII.   TRANSVERSE  POISSON'S  RATIO  STUDY 


Analysis 


Ef/Em 


V. 


V 


m 


V, 


1 
2 
3 
4 


1 

.7 

.25 

.25 

.25 

.250 

6 

.7 

.25 

.25 

.25 

.219 

120 

.7 

.25 

.25 

.25 

.109 

120 

.7 

.25 

.30 

.26 

.134 

The  results  of  this  study  indicate  that  the  Poisson's 
ratio  in  the  transverse  direction  is  also  greatly  dependent 
on  the  ratio  of  constituent  moduli  of  elasticity.   The 
greater  the  ratio  E  /E  ,  the  greater  the  variance  from  the 
Poisson's  ratio  of  the  constituent  materials. 

At  small  filament  volume  fractions  the  transverse 

Poisson's  ratio  exceeds  that  of  either  constituent.   At 

large  values  of  filament  volume  content  vT  is  considerably 

less  than  either  constituent.   At  50%  filament  volume  the 

transverse  Poisson's  ratio  is  approximately  equal  to  the 

average  of  the  constituent  Poisson's  ratios,  irrespective 

of  the  value  of  E^/E  .   This  is  shown  in  Fig.  (14). 

f   m 

Comparison  of  analyses  3  and  4  of  Table  VII  shows 

that  an  increase  in  constituent  Poisson's  ratio  does  effect 

vT  but  not  as  significantly  as  does  the  ratio  Ef/E  .   A  20% 

increase  in  v   resulted  in  approximately  a  20%  increase  in 

v  .   The  effect  of  changing  constituents  Poisson's  ratio  is 

limited  however  by  the  range  of  values  for  available 

materials . 
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m 


Figure  (14) .   Effect  of  Constituent  Stiffness  on  Transverse 
Poisson's  Ratio. 

In  the  course  of  attempting  to  correlate  the  results 

of  this  study  with  other  theoretical  results  several 

alternative  analytical  methods  were  considered.   Some  of 

these  were  found  to  give  results  for  Poisson's  ratio 

comparable  to  other  extreme  values  found  in  the  literature. 

These  alternative  methods  were  discarded  as  erroneous.   To 

assist  in  evaluating  the  various  predicted  values  of  vT 

two  of  these  methods  will  be  discussed  here. 

i)   If  it  is  assumed  that  the  composite  exhibits 

the  same  stiffness  in  all  directions  the  stress-strain 

relations  become 


£  =  -  (a  -  v_a  ) 

x     e     *       T  y 


=  i  (aw  -  Vm<0 


E 


T  x' 


L  1 

1   (ax  +  ay)  +  e  °z 


(111-22) 


Using  the  data  of  the  various  load  cases  these  equations 
yield  values  for  longitudinal  Poisson's  ratio,  v  ,  that 
follow  a  law  of  mixtures.   As  noted  in  the  previous  section 
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this  concurs  with  most  theoretical  values.  The  transverse 
Poisson's  ratio,  however,  takes  on  extremely  small  values, 
considerably  less  than  .02  for  60%  models.  This  method  is 
not  correct  due  to  the  inaccurate  assumption  regarding  the 
stress-strain  relations  of  the  composite. 

ii)   It  was  also  attempted  to  evaluate  the  transverse 
Poisson's  ratio  by  simulating  a  tension  test.   A  finite 
element  analysis  was  carried  out  on  the  model  shown  in  Fig. 
(15a) .   The  boundary  conditions  were 

u   =0   on   y  =  0 

y        * 

u  =0   on   x  =  0 

u  =  1   on  x  =  a 
This  resulted  in  a  displacement  along  the  boundary  as  shown 
in  Fig.  (15b) . 


— *» 

I""-1 

^. 

>J 

1 

J 

(a)   Model 


0    (b)   Boundary  Displacement 


Figure  (15).   Simulated  Tension  Test  of  a  Composite. 

Consequently  the  Poisson's  ratio,  v  ,  varys  along  the 
boundary.   In  the  problem  considered  the  values  ranged 
from  approximately  0.15  to  0.7.   The  error  in  this  method 
is  that  the  rectangular  boundaries  of  the  basic  block  do 
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not  remain  rectangular,  leading  to  violation  of 
compatibility  requirements. 

The  above  result  does  however  suggest  a  physical 
understanding  of  the  Poisson's  ratio  found  in  this 
investigation.   Figure  (16a)  demonstrates  how  the  results 
of  this  study  require  the  composite  to  act  under  a  tensile 
load.   Figure  (16b)  shows  how  the  material  would  react 
if  compatible  boundary  conditions  were  not  enforced.   In 
Fig.  (16c)  the  stresses  necessary  to  ensure  the  required 
compatible  boundary  conditions  are  shown.   The  stresses 
are  such  as  to  cause  the  displacements  curve  L?  to 
coincide  with  l-i  •   The  integral  of  these  stresses  is  zero, 


i.e. ,  S   =  0. 

y 

a  ^  (on  boundary) 


a  =0  (on  boundary)   ay  7*0  (on  boundary) 


f£> 


(b) 


-He  K 


(c) 


-^e  f< 


Figure  (16).   Transverse  Poisson's  Ratio. 

6 .   Stress  Concentration  Factor  (SCF) 

Before  examining  the  results  obtained  for  stress 
concentration  factors  some  consideration  must  first  be 
given  to  the  value  and  interpretation  of  such  results. 
The  primary  value  of  SCF  is  in  its  relationship  to  the 
strength  of  a  material.   For  a  homogeneous  material 
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initial  yield  will  occur  at  the  point  of  highest  micro- 
stress,  that  is  the  location  of  the  SCF.   For  identical 
materials  with  different  SCF,  which  may  occur  as  a  result 
of  loading  conditions,  initial  yield  will  occur  at  a  lower 
stress  level  in  the  material  with  the  larger  stress  con- 
centration factor. 

For  composite  materials  SCF  requires  further  interpre- 
tation.  Due  to  material  discontinuity,  stress  distributions 
will  differ  in  the  matrix  and  filament  materials.   For  this 
reason  SCF  should  be  calculated  separately  for  filament  and 
matrix  regions  of  a  composite.   If  the  matrix  is  the  larger 
no  problem  is  posed,  yield  will  occur  in  the  matrix  material. 
If  the  filament  values  are  larger  further  investigation  is 
necessary  to  ensure  that  initial  yield  is  not  occurring  in 
the  filament.   If  the  ratio  of  filament  to  matrix  stress 
concentration  factor  is  less  than  the  ratio  of  filament  to 
matrix  yield  strength,  yield  will  still  occur  in  the  matrix. 
That  is 


yield  in  matrix 


SCFf 

/v 

SCF  i 
m/ 

m 

<  1 

SCFf 
m/ 

m 

=  1 

SCFf 

/s 

>  1 

SCF  / 
m/ 

Yo 

m 

simultaneous  yield  in 
matrix  and  filament 


yield  in  filament     (111-23) 
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Caution  should  be  exercised  in  interpreting  stress 
concentration  factors.   Assuming  that  yield  occurs  in 
the  matrix  the  macroscopic  load  required  to  cause  yield 

may  be  related  to  the  SCF  as 

Y 

°m 

SCF  - (111-24) 

sc 

where  S   is  the  composite  yield  strength.   This  may  represent 

a  longitudinal,  transverse,  or  longitudinal  shear  yield 

strength.   For  different  materials  with  identical  matrix 

yield  strength,  YQ  ,  therefore  a  smaller  SCF  will  require 

m 
a  larger  macroscopic  load,  S  ,  to  cause  yield.   Thus  for 

such  materials  the  SCF  may  be  directly  related  to  yield 
strength.   An  appealing  aspect  of  this  result  is  that  for 
composites  having  identical  matrix  yield  values,  SCF  offers 
an  immediate  evaluation  of  comparative  yield  strengths  for 
the  longitudinal,  transverse,  and  longitudinal  shear  loading 
cases. 

For  any  composite  with  initial  yield  occurring  in  the 
matrix,  regardless  of  the  matrix  yield  strength,  the  compos- 
ite yield  strength  may  be  calculated  for  the  three  loading 
cases  as 

S   =  Y   /SCF  (111-25) 

c     °m 

m 

As  noted  above  SCF  also  allows  a  rapid  check  on  which 
constituent  first  reaches  yield. 

For  transverse  and  longitudinal  shear  loading  the 
maximum  SCF  usually  occurs  in  the  matrix  material.   For 
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this  reason  any  reference  to  SCF  for  these  specific  cases 
should  be  assumed  to  be  a  matrix  property  unless  specifically 
stated  to  the  contrary.   No  such  generalization  is  warranted 
for  the  longitudinal  loading  problem. 

Several  authors  have  treated  stress  concentration  factors 
extensively.   Adams  and  Doner ,  [3]  [6]    using  a  finite  difference 
technique,  have  provided  a  great  deal  of  numerical  data  for 
both  the  transverse  and  longitudinal  shear  loading  problems. 
For  the  transverse  problem  with  a  55%  Glass/Epoxy  composite 
Adams  and  Doner  calculated  SCFX  =  1.86.   They  also  presented 
comparative  data  from  solutions  utilizing  other  analytical 
methods,  ranging  in  value  from  1.49  to  1.78.   Analyses  were 
not  performed  here  for  55%  volume  composites,  however  a  47% 
model  gave  SCFX  =  1.49,  and  a  6  0%  model  gave  SCFX  =  1.70. 
This  data  is  included  in  Table  VIII. 

TABLE  VIII.   SCF  FOR  TRANSVERSE  NORMAL  LOADING 

E./E  =20     vf=.3   v  =.34   v.=.55 
f   m        £       m       f 

Method  of  Solution  Stress  Concentration  Factor,  SCFX 

Adams  and  Doner  [_6]  1.86 

Wilson     [l3]  1.71 

Pickett     |4]  1.49 

Foye    [7]  1.78 

Present  analysis  (47%)  1.49 

Present  analysis  (60%)  1.70 

Present  analvsis  (50%  estimate)  1.63 
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For  the  longitudinal  shear  problem  Adams  and  Doner  [3] 
gave  results  for  SCFX  of  40%  and  55%  models  as  1.605  and 
1.787  for  Glass/Epoxy  and  1.649  and  1.854  for  Boron/Epoxy. 
As  shown  in  Table  IX  these  values  bracket  the  present  47% 
solutions  of  1.66  for  Glass/Epoxy  and  1.69  for  Boron/Epoxy 

TABLE  IX.   SCF  FOR  LONGITUDINAL  SHEAR  LOADING 


Material 


Filament  Volume 


SCFX 


Source 


Glass/Epoxy 
Glass/Epoxy 
Glass/Epoxy 
Boron/Epoxy 
Boron/Epoxy 
Boron/Epoxy 


40% 
47% 
55% 
40% 
47% 
55% 


1.6  05  Adams  and  Doner 

1.66  Present  Analysis 

1.787  Adams  and  Doner 

1.649  Adams  and  Doner 

1.69  Present  Analysis 

1.854  Adams  and  Doner 
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IV.   PARAMETRIC  ANALYSIS  OF  COMPOSITES 

The  purpose  of  this  section  is  to  analyze  the  effect 
that  individual  constituent  material  properties  have  on 
the  mechanical  properties  of  the  composite.   Properties 
considered  include  the  modulus  of  elasticity,  E,  Poisson's 
ratio,  v,  and  yield  strength,  Y  ,  of  both  matrix  and 
filament  materials.   In  this  study  it  is  assumed  that  both 
filament  and  matrix  materials  are  isotropic. 

The  goal  is  not  to  investigate  every  possible  combination 
of  material  properties  but  rather  to  establish  criteria  and 
guidance  for  predicting  composite  properties  associated  with 
specific  constituent  properties. 

The  macroscopic  properties  that  are  considered  as  per- 
formance indices  include  the  initial  yield  strength,  the 
specific  strength,  longitudinal  and  transverse  stiffness, 
shear  modulus,  stress  concentration  factors  and  specific 
moduli. 

To  determine  how  constituent  properties  effect  composite 
behavior,  parameter  variations  were  carried  out  on  ficti- 
tious materials.   Variations  were  performed  by  starting  with 
some  actual  properties,  primarily  those  of  a  Boron/Aluminum 
or  Graphite/Epoxy  composite,  and  independently  varying  the 
particular  property  in  question. 
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A.   EFFECT  OF  FILAMENT  VOLUME 

To  evaluate  the  effect  of  filament  volume  on  composite 
properties  several  groups  of  40%,  50%,  and  60%  filament 
volume  models  were  investigated. .  These  included 

1)  Boron/NARMCO  Epoxy 

2)  High  modulus  graphite  (MODMOR  I)/4617  Epoxy 

3)  High  strength  graphite  (MODMOR  II)/NARMCO  Epoxy 

4)  MODMOR  11/4  617  Epoxy 

For  MODMOR  I/NARMCO  Epoxy  additional  models  of  3  0%  and  7  0% 
were  also  utilized.  The  material  properties  for  all  mate- 
rials are  listed  in  Table  II. 

In  a  composite  material  subjected  to  a  constant  longitu- 
dinal strain  the  macroscopic  stress  is  related  to  the  micro- 
scopic axial  stresses  as 

S   =  vfof  +  v„ct  (IV-1) 

c    £  f    m  m 

where  v_  and  v   represent  the  volume  fraction  of  filament 

f      m   c 

and  matrix  materials.   An  assumption  of  this  study  is  that 
the  filament  and  matrix  are  perfectly  bonded.   For  longitu- 
dinal loading  the  strain  will  be  constant  over  any  transverse 
cross  section,  thus; 

ec  -  H   =  Sn  (IV"2) 

Within  the  elastic  region  the  stress  strain  relations  give 

ET  =  v-rE.  +  v  E  (IV-3) 

L     f  f     m  m 

The  longitudinal  stiffness  of  a  composite  should  be  described 
then  by  a  linear  relation,  the  so  called  law  of  mixtures. 
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This  line  passes  through  the  matrix  modulus  at  0%  filament 

volume  (ET  =  E   when  v^  =  0)  and  through  the  filament 
L    m       f  3 

modulus  at  100%  filament  volume  (E   =  E   when  vf  =  1) ,  as 

Li        T.  E 

shown  in  Fig.  (17)  for  MODMOR  I/NARMCO  Epoxy. 

For  the  case  of  transverse  stiffness  the  relationship 
is  not  linear  with  respect  to  filament  volume.   For  loading 
in  the  X  or  Y  direction  a  constant  strain  applied  at  the 
boundary  will  not  result  in  a  uniform  strain  throughout  the 
material.   This  is  due  to  the  change  in  material  composition 
for  cross  sections  perpendicular  to  the  direction  of  loading. 
At  zero  filament  volume  the  material  has  the  properties  of 
the  matrix,  that  is  the  composite  transverse  modulus,  E  , 
must  equal  the  matrix  modulus,  E  .   At  100%  filament  volume 
the  material  will  have  the  properties  of  the  filament  and 
the  composite  transverse  modulus,  E  ,  will  equal  the  filament 
modulus,  Ef.   The  transverse  stiffness,  Em/  versus  filament 
volume,  vf,  curve  must  go  through 

ET  "  Em   at   vf  =  ° 
and  increase  monotonically ,  terminating  at 

Em  =  Ef   at  vf  =  1. 

All  cases  investigated  followed  this  pattern.   The  results 

for  MODMOR  I/NARMCO  Epoxy  are  plotted  in  Fig.  (17).   The 

increasing  slope  of  the  curve  indicates  that  an  increase  in 

filament  volume  has  a  more  significant  effect  on  transverse 

stiffness  at  higher  filament  volume  percentages.   For  example, 

for  a  MODMOR  I/NARMCO  Epoxy  composite  a  15%  increase  in 

filament  volume,  v^,  from  35%  to  50%  results  in  a  42%  increase 

in  composite  transverse  stiffness,  ET.   A  similar  15%  increase 
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Figure  (17).   Composite  Elastic  Properties  vs.  Filament  Volume 
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in  v   from  60%  to  75%  results  in  a  75%  increase  in  composite 
transverse  stiffness.   For  composite  shear  modulus,  G,  a 
change  in  cross  sectional  material  properties  indicates  a 
behavior  similar  to  that  of  the  transverse  stiffness. 

In  longitudinal  loading  problem  III,  the  stress  con- 
centration factor  SCFZ  is  calculated.   The  unit  strain 
applied  results  in  a  uniform  stress  in  the  matrix  material 
and  a  uniform  stress  of  a  different  magnitude  in  the 
filament  material.   For  the  same  constituent  material  and 
identical  unit  strain  these  stresses  should  not  vary  with 
changes  in  filament  volume  since  they  must  satisfy  the 
relations 

of   =   Efe  (IV-4) 

am  =  V  <IV"5> 

The  resultant  surface  traction,  S  ,  the  average  value  of 

the  microstresses  o ^    and  a   acting  on  the  cross  section, 

£      m 

will  change  proportionally  with  the  increase  in  filament 
volume,  according  to  the  equation 

S   =  vPo.   +   v  om  (IV-6) 

z     f  f     m  m  v     ' 

Thus  as  (az)    ,  the  numerator  of  SCFZ,  Eq.  (III-ll) , 
max 

remains  constant,  the  denominator  S   increases  lineally 
with  increasing  filament  volume,  v f ,  and  the  stress  con- 
centration factor  itself  decreases. 

For  a  homogeneous  material  the  uniform  stress  throughout 
the  entire  block  would  equal  the  applied  stress  and  SCFZ 
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would  equal  one.   The  curve  for  SCFZ  versus  filament 
volume  should  therefore  pass  through  the  points 
SCFZ  (vf  =0)  =1,  and  SCFZ  (vf  =1)  =1. 

Consider  the  case  of  the  SCFZ  in  the  filament.   The 
addition  of  a  small  amount  of  filament  material  to  a  0% 
model  would  result  in  a  negligible  change  in  the  surface 
traction.   That  is,  for  unit  strain 


e   =  1.0  (IV-7) 

z 

S   ~  a    =  Eme  =  Em  (IV-8) 

Z   m    m     m 

The  maximum  filament  stress  in  this  case  would  be 


oz      =  Efe  =  Ef  (IV-9) 

For  a  small  filament  volume  the  filament  SCFZ  would  be 

°f        Ef 
SCFZ  =  —  =  _JL  (IV-10) 

Om         E 
m    m 

Figure  (18)  for  Boron/6061  Aluminum  shows  this  graphically. 

The  points  at  0%  and  100%  filament  volume  are  established  in 

accordance  with  the  above  discussion  while  the  other  points 

are  results  of  the  analysis. 

For  the  matrix  SCFZ  the  maximum  stress  will  remain 

constant  while  the  average  surface  traction  increases  with 

increase  of  vf  ,  as  noted  in  Eq.  (IV-6) .   For  the  matrix  case 

therefore  the  curve  will  originate  at  SCFZ  (v   =0)  =1, 

and  decrease  to  the  point  where,  with  little  matrix  material 

remaining, 

E 
SCFZ  =  =21  (IV-11) 

bf 
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Recall  from  section  II-B-6  that  if  the  ratio  of  the 

filament  SCF  to  the  matrix  SCF  is  less  than  the  ratio  of 

filament  to  matrix  yield  strengths,  YQ  /Y   ,  yield  will 

f    m 
occur  in  the  matrix.   This  will  determine  in  each  individual 

case  which  of  the  stress  concentration  factors  discussed 

above  is  pertinent. 

SCFX  and  SCFXZ  increase  directly  with  filament  volume, 

as  shown  in  Fig.  (18).   These  stress  concentrations  occur 

in  the  matrix  at  the  interface. 

B.   EFFECT  OF  FILAMENT  YOUNG'S  MODULUS 

To  evaluate  the  effect  of  variations  in  the  modulus  of 
elasticity  of  the  filament,  Ef,  on  the  composite  mechanical 
properties  four  computer  analyses  were  made.   Each  analysis 
used  a  50%  filament  volume  model  with  material  properties 
as  listed  in  Table  X. 

TABLE  X.   FILAMENT  YOUNG'S  MODULUS  VARIATION 

Ef  °f  Em  °m 

Analysis   (psixlO"6)   vf    (psixl0~5)   (psixlO-6)   v   (psixlO-1*) 

1  30  .2  4.8  10.  .3  .3 

2  40  .2  4.8  10.  .3  .3 

3  60  .2  4.8  10.  .3  .3 

4  70  .2  4.8  10.  .3  .3 
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Figure  (19)  displays  the  transverse  stiffness,  longitu- 
dinal stiffness,  and  shear  modulus  of  the  composite  as  a 
function  of  fiber  modulus.   This  figure  shows  that  the  most 
significant  effect  of  increasing  the  filament  modulus  is 
the  linear  increase  in  the  composite  longitudinal  stiffness, 
E  ,   This  relates  with  matrix  modulus  variation  and  will  be 

Xj 

treated  in  greater  detail  in  the  following  two  sections. 

The  relationship  between  the  composite  longitudinal 

stiffness,  E  ,  and  the  filament  Young's  modulus,  E^,  is 
L  t 

linear.   In  fact  for  the  same  reasons  as  discussed  in 
section  IV-A  for  filament  volume  variation,  the  relation- 
ship can  be  accurately  described  by  a  rule  of  mixtures  as 

E   =  v^E.  +  v  E  (IV-12) 

L     f  f     m  m 

The  plot  of  the  composite  shear  modulus  variation  of 
Fig.  (19)  appears  to  be  linear,  although  it  will  be  shown 
that  this  variation  must  be  non-linear.   Consider  that,  for 
the  combination  of  materials  considered  here,  when  the 
composite  is  homogeneous  E   =  10  x  106psi.   Hence 

G   a  Q  s  g   =  3.84  x  106  psi 
c    f    m  r 

The  shear  modulus  curve  of  Fig.  (19)  must  vary  considerably 
from  its  apparent  linear  form  to  pass  through  this  point, 
as  shown  in  Fig.  (25). 

As  would  be  expected,  considering  the  linear  variation 
of  longitudinal  stiffness,  the  longitudinal  yield  strength 
also  increased  linearly  with  filament  modulus  increase. 
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Figure  (19).   Composite  Elastic  Properties  vs. 
Filament  Modulus  of  Elasticity. 
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Variation  in  transverse  and  longitudinal  shear  strength 
was  negligible  for  the  range  of  E   considered,  E   =  30-75  x  10epSi 

C.   EFFECT  OF  MATRIX  YOUNG ' S  MODULUS 

To  investigate  the  effect  of  variation  of  the  matrix 
modulus  of  elasticity  on  the  composite  mechanical  properties 
a  series  of  four  analyses  were  made.   All  analyses  were  made 
with  a  50%  filament  volume  model  with  constituent  properties 
as  shown  in  Table  XI . 

TABLE  XI.   MATRIX  YOUNG'S  MODULUS  VARIATION 


Ana! 


m 


m 


-  5 


1-6 


lysis   (psixlO-6)   vf   (psixlO  5)   (psixlO  6)  vm   (psixlO  5) 


1 

60. 

.2 

4.8 

2 

60. 

.2 

4.8 

3 

60. 

.2 

4.8 

4 

60. 

.2 

4.8 

10.  .3  .346 

7.  .3  .346 

5.  .3  .346 

2.5  .3  .346 


As  established  in  Eq.  (IV-3)  the  composite  longitudinal 
stiffness  varys  lineally  as 


ET  =  v-rE^  +  v  E 
Jj     r  f     mm 


Thus  the  variation  of  E   with  matrix  modulus  is  also  linear 

L 

as  shown  in  Fig.  (20) .   From  this  figure  it  is  seen  that  the 

variation  of  E   has  a  greater  effect  on  the  transverse 
m 

stiffness  and  shear  modulus  than  on  longitudinal  stiffness. 
An  increase  of  400%  in  the  matrix  modulus  of  elasticity 
results  in  percentage  increases  of 
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Figure  (20).   Composite  Elastic  Properties  vs. 
Matrix  Modulus  of  Elasticity. 
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i)   E   =  14% 

Xj 

ii)   ET  =  207% 
iii)    G  =  222% 
This  also  resulted  in  an  increase  in  transverse  and  longitu- 
dinal shear  yield  strength.   The  complementary  effect  of 
this  result  with  that  of  the  filament  modulus  variation 
will  be  discussed  in  the  following  section. 

An  interesting  aspect  of  the  variation  of  the  matrix 
Young's  modulus  is  that,  despite  the  increase  in  longitudinal 
stiffness,  the  longitudinal  yield  strength,  SL,  decreases. 
The  reason  for  this  may  be  seen  in  Fig.  (21). 
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Figure  (21).   Effect  of  Matrix  Stiffness  on  Longitudinal 
Yield  Strength. 

For  two  matrix  materials  with  similar  initial  yield  points 
one  with  the  lower  Young's  modulus  will  require  a  larger 
\   strain  to  reach  yield.   Due  to  the  relatively  small  effect 
of  matrix  stiffness  on  composite  longitudinal  stiffness,  E  , 
assume  that  E   is  the  same  for  both  cases.   The  larger  strain, 

Li 
£2/  will  result  in  a  larger  composite  yield  strength,  S-j-  . 
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D.   EFFECT  OF  CONSTITUENT  PROPERTY  RATIO 

The  results  noted  for  filament  and  matrix  stiffness 
variation  may  be  conveniently  correlated  through  use  of 
non-dimensional  parameters.   For  example,  composite  longi- 
tudinal stiffness,  E  ,  may  be  normalized  with  respect  to 

JL 

filament  stiffness,  i.e.,  E  /Ef,  and  plotted  against  the 

constituent  modulus  ratio  E_/E^.   Figure  (22)  is  such  a 

m   r      3 

representation.   Longitudinal  stiffness  may  be  expressed  as 

EL  =  vfEf  +  vmEm  (IV-13) 

Dividing  through  by  Ef  and  rearranging  terms  yields 

EL/Ef  =  vf  (1  -  Em/Ef)  +  Em/Ef  (IV-14) 

Equation  (IV-14)  describes  the  configuration  of  the  non- 
dimensional  longitudinal  stiffness  curve.  The  following 
points  are  significant: 

i)   when  E  —  E~ ,    i.e.,  a  homogeneous  material, 

EL/Ef  =  VEf  "  l-  (IV"15) 

ii)   when  E   =  0,  i.e.,  zero  matrix  stiffness, 

EL/Ef  =  vf (1)  =  v  (IV-16) 

Equations  (IV-15)  and  (IV-16)  were  used  to  determine  the 
extreme  points  of  Fig.  (22) .   The  points  noted  on  this 
figure  are  the  results  of  this  investigation  for  a  filament 
volume  fraction,  v^  =  .5.   For  comparison  some  other  curves 
were  drawn  using  the  endpoints  determined  by  Eqs.  (IV-15)  and 
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(IV-16) .   The  non-dimensional  longitudinal  stiffness 
curve  for  any  filament  volume  must  fall  in  the  triangular 
region  of  Fig.  (22)  formed  by  the  Vf  =  0  curve  and  the 
vf  =  1.0  curve. 

The  comparative  effect  of  filament  Young's  modulus 
and  matrix  Young's  modulus  variations  on  composite  longi- 
tudinal stiffness  may  be  examined  by  consideration  of  Eq. 
(IV-13) .   Taking  partial  derivatives  of  this  equation  with 
respect  to  the  individual  independent  quantities  yields 


8EL 

=  v  (IV-17) 

3E      m 
m 

and 


dE 
L 

8Ef 


vf  (IV-18) 


Thus  the  effect  of  constituent  stiffness  variation  on  longi- 
tudinal stiffness  is  directly  related  to  the  volume  fraction 
of  the  constituent.   For  the  50%  filament  volume  case  this 
is  seen  to  be  true  as  Figs.  (19)  and  (20)  show  that  filament 
and  matrix  stiffness  variation,  respectively,  result  in 
identical  slopes  of  0.5  when  plotted  versus  longitudinal 
stiffness. 

These  comments  concerning  comparative  effects  of  matrix 
and  filament  stiffness  variation  refer  to  the  change  per  unit 
stiffness  of  constituent.  The  actual  range  of  longitudinal 
stiffness  due  to  matrix  and  filament  stiffness  variation  is 
considerably  different  however.  Consider  Fig.  (23),  a  non- 
dimensional  longitudinal  stiffness  curve  for  some  arbitrary 
filament  volume  fraction,  v^ . 
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Figure  (22) 


Effect  of  Constituent  Stiffness  on 
Longitudinal  Stiffness. 
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Figure  (23).   Non-Dimensional  Longitudinal  Stiffness. 

For  a  given  constituent  moduli  ratio,  say  point  A,  a 
specific  ratio,  E  /E  ,  results  at  point  A*.   At  this  point 

Li    I 

EL   =  A'E 
A        A 

If  the  constituent  ratio  is  increased  to  point  B  by  increas- 
ing the  matrix  modulus  while  holding  the  filament  modulus 
constant  then 

A      B 

ET   =  B'E. 
LB       rB 

The  ratio  of  longitudinal  stiffness  for  these  two  hypothetical 

composites  is 


B'E. 


B 


B 


(ET/E  ) 
B'      L   f  B 


A'E, 


A' 


A 


(EL/Ef) 


(IV-19) 


Thus  variation  in  the  matrix  modulus  results  in  a  change  in 
composite  longitudinal  stiffness  in  direct  proportion  to  the 
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change  in  non-dimensional  stiffness,  E  /Ef.   In  general 

Li    •*- 

the  maximum  range  of  non-dimensional  stiffness,  E  /E,.,  is 

J_i     L 

from  the  value  of  filament  volume  fraction,  vf ,  to  1.0. 

The  maximum  increase  in  longitudinal  stiffness  E  ,  due 

L 

solely  to  a  variation  of  matrix  modulus  is 

\ 


% 


EL 


1.0 

'  max 


For  the  50%  filament  volume  case  of  Fig.  (22)  this  maximum 
increase  is  a  factor  of 


1.0 

=  2.0 


.5 


If  the  constituent  ratio  of  Fig.  (23)  is  increased  to 
point  B  by  decreasing  the  filament  modulus  while  holding 
the  matrix  modulus  constant,  then 


Ef   <   Ef 
B       A 


ET   =  A'E.       ET   =  B'E_ 
LA       fA       LB       fB 


EL     A'Ef 

(IV-20) 


E      B'E 
LB       fB 


Thus  the  increase  in  longitudinal  stiffness  with  variation 
of  the  filament  modulus  is  dependent  on  the  non-dimensional 
stiffness,  E  /Ef,  and  the  ratio  of  filament  moduli,  Ef  /Ef 


B 
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For  actual  filament  materials  the  range  for  the  filament 
modulus  is  approximately  10  x  106  to  60  x  106  psi.   Consider 
a  constant  matrix  material,  for  example  aluminum.   For  the 
50%  filament  volume  case  of  Fig.  (22)  the  range  of  E  /Ef 
corresponding  to  this  range  of  filament  modulus  is  50  to 
1.0.   For  this  50%  case  the  variation  in  longitudinal  stiff- 
ness due  to  variation  in  ther^f i lament  modulus  alone  is 


\ 


EL 
B 

Variation  of  the  filament  modulus  is  thus  seen  to  result  in 

a  much  greater  range  of  longitudinal  stiffness  than  does 

the  variation  of  the  matrix  modulus  of  elasticity  because 

of  the  greater  range  of  E   as  compared  to  E  . 

Figure  (24)  is  a  non-dimensional  representation  of  the 

effect  of  constituent  stiffness,  Em/Ef ,  on  a  normalized 

composite  transverse  stiffness,  E_/Ef.   Note  that  the  curve 

approaches  E   =  0  as  Em/E   approaches  zero.   This  indicates 

that  when  the  matrix  material  is  extremely  soft  in  comparison 

to  the  filament  that  the  macroscopic  transverse  stress 

required  to  produce  a  unit  displacement  of  the  block  boundary, 

as  in  load  case  I,  is  virtually  zero.   In  the  limit,  when 

E   =0,  then  E   =0.   Thus  it  is  theoretically  possible  to 
m  T 

produce  a  composite  with  a  value  of  transverse  stiffness 
approaching  zero. 

Comparison  of  Figs.  (19)  and  (20)  shows  that  the  slope 

of  the  Em  vs.  matrix  stiffness  curve  is  greater  than  that 

T  * 

80 
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m 


Figure  (24)  . 


Effect  of  Constituent  Stiffness  on 
Transverse  Stiffness. 
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of  the  E  vs.  filament  stiffness  curve.   Therefore  on  a 
T 

unit  stiffness  basis,  variation  of  matrix  Young's  modulus 
is  more  effective  than  variation  of  filament  Young's  modulus 
on  the  composite  transverse  stiffness. 

By  considering  a  non-dimensional  transverse  stiffness 
curve  and  following  the  same  technique  as  for  the  non- 
dimensional  longitudinal  stiffness  the  range  of  the  trans- 
verse stiffness  due  to  matrix  Young's  modulus  variation, 
analogous  to  Eq.  (IV-19) ,  is 


B 


max 


(ET/Ef) 


B 


(ET/Ef) 

A 


(IV-21) 


max 


The  range  of  E  /Ef/  as  is  seen  in  Fig.  (24),  is  zero  to  one. 
Thus  the  theoretical  maximum  increase  of  transverse  stiff- 
ness due  to  matrix  Young's  increase  is 


O 

T 
B 


Er 


1.0 


max 


The  maximum  variation  in  transverse  modulus  due  to  variation 
of  the  filament  stiffness  alone  is 


(IV-22) 


ET   \ 

A 

max 

(E„/E    ) 

A 

E*A 

Ef 
B 

Jjrp 

Bi 

(ET/Ef) 

B 

max 


For  common  filament  materials  the  range  of  filament  modulus 
is  approximately  10  x  106  to  60  x  106  psi.   Consider  a 
constant  matrix  material,  for  example  aluminum.   For  the  50% 
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filament  volume  case  of  Fig.  (24)  the  range  of  E  /Ef 
corresponding  to  this  range  of  filament  modulus  is  .349 
to  1.0.   Thus  the  maximum  variation  in  transverse  stiff- 
ness due  to  variation  of  the  filament  modulus  alone  is 


ET         1.01. 
i  B7  max  \    /  \ 


10.. 


The  range  of  transverse  stiffness  due  to  matrix  Young's 
modulus  variation  is  greater  therefore  than  that  due  to 
filament  Young's  modulus  variation. 

Figure  (25)  is  a  plot  of  the  normalized  composite  shear 

modulus,  G/Gr:,    versus  the  ratio  of  constituent  shear  moduli, 

G  /Gc.      This  curve  substantiates  earlier  contentions  that 
m'  f 

shear  modulus  variation  with  respect  to  constituent  moduli 
is  non-linear.   Again  this  curve  approaches  zero  as  the 
constituent  moduli  ratio  approaches  zero. 

Again  comparing  Figs.  (19)  and  (20)  it  is  seen  that 
the  slope  of  the  composite  shear  modulus  versus  matrix 
stiffness  is  greater  than  the  slope  of  the  composite  shear 
modulus  versus  filament  stiffness  curve.   This  indicates 
that,  per  unit  of  constituent  stiffness,  the  matrix  material 
has  a  greater  effect  on  composite  shear  modulus. 

By  considering  a  non-dimensional  shear  modulus  curve 
and  following  the  same  technique  as  for  the  non-dimensional 

longitudinal  stiffness,  the  range  of  shear  modulus  variation 
due  to  matrix  shear  modulus  variation,  analogous  to  Eq.  (IV-19) , 
is 
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B 


\  A 


max 


(G/Gf) 


B 


(G/Gf) 


(IV-23) 


max 


The  range  of  G/Gf ,  as  seen  in  Fig.  (25),  is  from  zero  to 

1.0.   Thus  the  theoretical  maximum  value  of  Eq.  (IV-23)  is 

/ 
G_  \ 

1.0 


^B 


\ 


0 


A 


max 


If  the  ratio  G  /Gf  is  increased  by  changing  the  filament 
modulus  the  maximum  range  attainable  is 

/ 


\ 


B  /  max 


(G/Gf)   Gf 
K        A 


(G/Gf)   Gf 
B    B 


(IV-24) 


max 


For  common  filament  materials  the  range  of  shear  modulus  is 
approximately  4.0  x  106  to  25.0  x  106  psi.   For  a  constant 
matrix  material,  for  example  aluminum,  the  range  of  G/Gf  of 
Fig.  (25)  corresponding  to  this  range  of  filament  shear 
modulus  is  .31  to  1.0.   These  values  are  for  a  50%  filament 
volume  case.   For  this  case  the  Eq.  (IV-24)  becomes 

A 


\GbI 


max 
Variation  of  the  matrix  shear  modulus  thus  results  in 

a  greater  range  of  composite  shear  modulus  than  does 

variation  of  filament  shear  modulus. 

The  effect  of  Poisson's  ratio  will  be  discussed  in  the 

next  section. 
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Figure  (25) . 


Effect  of  Constituent  Shear  Modulus  on 
Composite  Shear  Modulus . 
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The  primary  advantage  of  the  non-dimensional  represen- 
tations is  that  valid  comparisons  may  be  made  between 
composites  without  necessarily  holding  the  Young's  modulus 
of  one  constituent  constant. 

It  should  be  noted  that  the  non-dimensional  plots 
presented  here  are  for  a  50%  filament  volume  model.   Dif- 
ferent filament  volume  models  will  result  in  different 
curves. 

E.   EFFECT  OF  POISSON'S  RATIO 

Two  series  of  computer  analyses  were  made  for  evaluation 
of  the  effects  on  composite  behavior  due  to  change  of  Poisson's 
ratio  of  the  filament,  v^.   Four  analyses  were  made  using  a 
fictitious  material  approximating  the  properties  of  Boron/ 
Aluminum  while  varying  the  filament  Poisson's  ratio  from  .15 
to  .30  in  increments  of  .05. 

Another  series  of  analyses  was  conducted  with  MO DMO R  I 
graphite  and  epoxy.   This  was  done  because  Poisson's  ratio 
for  graphite  filaments  is  not  accurately  specified  and  it 
was  desired  to  determine  whether  or  not  this  would  signif- 
icantly effect  the  results  of  this  investigation. 

In  both  cases  no  significant  effect  was  noted  on  any 
composite  mechanical  property.   All  observed  variations 
were  less  than  1%. 

The  effect  of  the  variation  of  matrix  Poisson's  ratio 
was  evaluated  through  a  series  of  three  analyses.   Constituent 
properties  similar  to  Boron/6  061  Aluminum  were  used  while 
varying  the  matrix  Poisson's  ratio  from  to  .30  to  .35.   Again 
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no  significant  effect  on  composite  mechanical  properties 
was  noted.   It  is  expected  however  that  changes  in  Poisson's 
ratio  will  produce  more  significant  changes  for  post  yield 
behavior. 

F.   EFFECT  OF  MATRIX  YIELD  STRENGTH 

To  investigate  the  effect  of  the  matrix  yield  strength, 

Y   ,  on  composite  mechanical  properties,  a  series  of  four 

m 
computer  analyses  were  made.   All  used  a  50%  Boron  filament 

volume  model  with  constituent  properties  as  listed  in  Table 

XII. 


TABLE  XII.   MATRIX  YIELD  STRENGTH  VARIATION 


Ana.' 


1 
2 
3 

4 


of 


m 


-  5 


>-6 


m 


lysis   (psixlO  b)   vf   (psixlO  s )   (psixl0~b)  vm  (psixl'O  5) 


60.  .2  4.8 

60.  .2  4.8 

60.  .2  4.8 

60.  .2  4.8 


10.  .3  .346 

10.  .3  .230 

10.  .3  .130 

10.  .3  .030 


The  yield  point  of  a  material  has  no  effect  on  the 
elastic  properties  of  that  material.   It  merely  defines  the 
limit  of  elastic  behavior.   Thus  the  variation  of  matrix 
yield  strength  had  no  effect  on  the  composite  stiffness  or 
shear  modulus.   This  is  as  expected  considering  the  compar- 
ison in  section  III-B-4  of  two  aluminum  alloy  matrix 
materials . 
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Figure  (26)  shows  hypothetical  stress-strain  curves  for 

the  composite  (curve  A)  and  the  matrix  (curve  B) .   For  a 

given  matrix  yield  strength,  say  Yq  ,  a  specific  composite 

m 
macroscopic  yield  strength,  si,,  is  defined.   If,  as  stated, 

the  yield  point  of  a  material  does  not  effect  the  elastic 

properties,  then  the  stress-strain  curves  remain  the  same, 

independent  of "the  yield  strength  variation.   For  a  given 

increase  in  matrix  yield  strength,  say  a  50%  increase  to 

Y*  ,  plane  geometry  shows  that  an  identical  increase,  on 

m 
a  percentage  basis,  must  occur  in  the  composite  yield 

strength.   Figure  (27)  confirms  this  hypothesis. 


Figure  (26) .   Effect  of  Matrix  Yield  Strength  Variation, 

It  should  be  noted  that  the  above  discussion  assumes 
that  yield  occurs  in  the  matrix,  as  is  true  for  most,  but 
not  all,  composites. 

The  implication  of  this  phenomenon  is  that  there  is  no 
variation  in  the  shape  of  the  initial  yield  surface  due 
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Figure  (27) .   Effect  of  Matrix  Yield  Strength  on 
Composite  Yield  Strength. 
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solely  to  variation  in  the  matrix  yield  strength.   That  is, 
the  initial  yield  surface  changes  in  size  only.   This  can 
be  seen  in  Figs.  (C-l)  and  (C-2)  of  Appendix  C. 

G.   EFFECT  OF  FILAMENT  CROSS  SECTION 

The  majority  of  the  computations  included  in  this  investi- 
gation are  based  on  longitudinal  filaments  with  a  circular 
cross  section.   This  was  done  because  it  is  assumed  that  this 
particular  cross  section  would  be  the  easiest  to  produce.   The 
question  arises  as  to  the  effect  change  of  shape  of  filament 
cross  section  would  have  on  composite  properties.   Such 
changes  might  be  the  result  of  fabrication. 

To  provide  an  answer  to  this  question  two  elliptical 
cross  sectional  filaments  were  investigated.   These  were 
identical  ellipses  with  major/minor  axes  ratio  of  1.25  to 
1.0.   One  model  had  the  major  axis  oriented  along  the  x-axis 
and  one  along  the  y-axis,  as  seen  in  Fig.  (28). 


x  axis      Major  Axis        y  axis 
Figure  (2S).   Elliptical  Cross  Sections. 


The  rule  of  mixtures  and  the  comments  of  section  IV-A 
leading  to  its  use  are  independent  of  cross  section  shape. 
Elliptical  and  circular  cross  section  filaments  of  the  same 
volume  should  therefore  produce  identical  results  for  the 
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longitudinal  composite  stiffness,  E  .   In  the  three  cases 

investigated,  50%  models  of  Boron/6061  Aluminum,  Boron/NARMCO 

Epoxy,  and  S-glass/NARMCO  Epoxy,  this  was  found  to  be  true. 

The  variation  in  E  was  less  than  1%  in  all  cases,  as  may  be 

Li 

seen  from  Table  XIII. 


TABLE  XIII.   EFFECT  OF  -FILAMENT  CROSS  SECTION 
ON  LONGITUDINAL  STRENGTH  AND  STIFFNESS 


Circular 

E_(psixlO*6)  S  (psixlO-4) 
L  L 


Elliptical 

E  (psixlO-6)  S  (psixlb"1*) 
L  L 


Boron/6061  Aluminum    37.08       11.4 
Boron/NARMCO  Epoxy     28.23       16.71 
S-glass/NARMCO  Epoxy    6.05        3.47 


37.25  10.1 

28.57  14.5 

6.105  3.04 


The  transverse  composite  stiffness  reacted  as  expected. 
In  all  cases  E   increased  in  the  direction  of  the  major  el- 
liptical axis  and  decreased  in  the  direction  of  the  minor 
axis  as  noted  in  Table  XIV.   This  seems  reasonable  since  more 
of  the  stiff  filament  material  is  oriented  in  the  direction 
of  the  major  axis.   On  a  percentage  basis  the  increase  in 
transverse  stiffness  in  the  direction  of  the  major  axis  was 
approximately  twice  that  of  the  decrease  in  the  direction  of 
the  minor  axis. 

An  interesting  result  of  the  cross  section  variation 
that  is  evident  in  Table  XIV  is  that  the  transverse  strength 
changed  in  a  manner  opposite  to  that  of  the  transverse  stiff- 
ness.  That  is,  the  transverse  strength  increased  in  the 
direction  of  the  minor  axis  and  decreased  in  the  direction 
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of  the  major  axis.   This  result  may  be  understood  by 
considering  the  micromechanics  of  the  problems  in  question. 
For  ease  of  discussion  call  the  transverse  load  case  with 
the  major  axis  oriented  along  the  x-axis  load  case  IA,  and 
that  with  the  major  axis  along  the  y-axis  as  load  case  IB. 
These  load  cases  are  shown  in  Fig.  (29).   Using  this 
notation,  the  -analysis  results  require,  for  transverse 
normal  loading,  at  yield 


sx  <   sx 


(IV-25) 


A      B 

The  definition  of  macrostresses  requires  that  the  integral 
of  the  microstresses  along  the  boundary  be  related  as 


k< 


B 


(IV-26) 


TABLE  XIV.   EFFECT  OF  FILAMENT  CROSS  SECTION 
ON  TRANSVERSE  STRENGTH  AND  STIFFNESS 


Circular 


Elliptical 
Major  Axis     Minor  Axis 


Boron/6061  Aluminum 
Boron/NARMCO  Epoxy 
S-glass/NARMCO  Epoxy 


ET  (psixlCT6) 


24.4    3.0 
1.64    .24 
1.52    ,245| 

ST  (psixlti"1*) 


26.6   3.5 
2.05   .19 
1.83   .195 


23.0   2.5 
1.46   .30 
1.38   .295 


Figure  (29)  shows  the  stress  and  strain  distribution  on 
the  x-boundary  (x  =  a)  for  a  unit  displacement  for  load 
cases  IA  and  IB.   In  both  cases  yield  occurs  in  the  matrix 
at  a  point  near  the  interface,  along  the  x-axis.   Both  cases 
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also  show  little  variation  in  stress  or  strain  along  the 
x-axis  between  the  interface  and  the  boundary.   For  example, 
for  load  case  A,  i.e.,  major  axis  along  the  x-axis,  problem 
I  for  a  50%  S-glass/Epoxy  model  gives  a  stress,  a,  at  the 
interface  of  .664  x  106  psi.   At  the  boundary  (x  =  a)  ax 
equals  .688  x  106  psi.   The  same  material  for  load  case  B 
gave  an  interface  value  of  <rv  as  .366  x  106  psi  and  a 
boundary  value  of  .389  x  106  psi.   Thus  the  strain  to  cause 
yield  will  be  approximately  the  same  at  the  yield  point,  i.e., 
the  interface,  and  at  the  boundary  (x  =  a,  y  =  0) .   If  the 
strain  and  stress  curves  of  Fig.  (29)  are  scaled  to  give 
identical  strain  values  at  the  boundary  (x  =  a,y  =  0)  the 
more  uniform  distribution  of  case  IB  will  result  in  a  larger 
value  of  the  integral  of  the  stress.   This  confirms  Eqs. 
(IV-25)  and(IV-26). 

It  seems  therefore  that  the  stiffness  of  the  composite 
depends  on  the  degree  of  filament  material  alignment  in  a 
particular  direction  while  the  strength  depends  on  the 
maximum  filament  cross  section  in  a  direction  perpendicular 
to  that  of  the  load.   The  inference  of  this  is  that  the 
optimum  cross-section  shape  would  approximate  that  of  a 
plane  whose  orientation  depends  on  the  desired  properties. 
This  is  shown  in  Fig.  (30) .   The  filament  shape  can  not  be 
a  rectangle  as  the  square  ends  would  induce  prohibitive 
stress  concentrations. 
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Figure  (3D) .   Approximate  Optimum  Cross  Section  for 
Particular  Properties. 


In  the  case  of  circular  cross  section  filaments  the 
transverse  properties  are  the  same  in  the  y-direction  as 
in  the  x-direction.   Thus  transverse  yield  strength  in 
the  y-direction  would  be  numerically  equivalent  to  that  in 
the  x-direction.   This  is  shown  in  Fig.  (31a).   Again  assume 
the  convenient  designation  of  load  case  A  as  that  of  the 
elliptical  cross  section  with  the  major  axis  along  the 
x-axis,  and  load  case  B  with  the  major  axis  along  the  y- 
axis.   The  results  of  these  two  load  cases  may  be  interpreted 
as  transverse  strength  in  different  directions  of  the  same 
model  as  shown  in  Fig.  (31b) . 

On  a  percentage  basis  the  change  in  transverse  strength 
due  to  cross  section  variation  was  considerable,  ranging  from 
15%  for  Boron/Aluminum  to  23%  for  Boron/NARMCO  Epoxy.   Of 
more  significance  however  is  the  fact  that  the  vector  sum 
of  Sv  and  S   is  independent  of  the  filament  cross  section. 

X        y  c 

This  indicates,  in  vector  notation,  for  Fig.  (31) 


S  -  Sx   +  Sy  "  Sx   +  Sy 
c     c     e     e 


(IV-27) 
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Or,  alternatively, 


(IV-28) 


<l       Circular  Filament 
~    (a) 


<=F> 


Ha  K- 


■Mah" 


t 


S    ?*  S„ 

e     e 

Elliptical  Filament,  Identical  vF 
(b) 

Figure  (31) .   Transverse  Strength  Relation  to  Filament 
Cro;G  Section. 

This  indicates  that  any  increase  in  strength  in  one  direction 
is  directly  offset  by  a  decrease  of  the  strength  in  the  per- 
pendicular direction. 

In  each  case  there  was  a  noticeable  increase  in  the  stress 
concentration  factors  SCFX  and  SCFXZ ,  although  SCFZ  remained 
relatively  constant.   The  net  result  was  a  reduction  in  the 
longitudinal  yield  strength  of  the  elliptical  cross  section 
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compared  to  the  circular  cross  section.   The  decrease  was 
approximately  12%  in  each  case.   The  implication  is  that  a 
circular  cross  section  is  the  most  effective,  with  respect 
to  a  given  filament  volume,  in  terms  of  longitudinal  strength, 
Adams  and  Doner  [6  J  reported  that  stress  concentration  factors 
are  directly  related  to  the  maximum  radius  of  curvature  of 
the  filament.  -For  a  given  filament  volume  a  circular  cross 
section  results  in  the  smallest  possible  maximum  radius  of 
curvature.   This  also  implies  that  the  circular  cross  section 
is  the  most  efficient. 

Another  significant  result  evident  in  the  analysis  of 
the  elliptical  filament  is  the  effect  of  filament  shape  on 
the  composite  shear  modulus,  G.   Normalized  with  respect  to 
the  circular  composite  shear  modulus  the  results  obtained 
are  listed  in  Table  XV. 


TABLE  XV.   EFFECT  OF  CROSS  SECTION  VARIATION 
ON  COMPOSITE  SHEAR  MODULUS 


Gell/Gc         f 
Composite  Major  Axis  Orientation 


Boron/6061  Aluminum 
Boron/NARMCO  Epoxy 
S-glass/NARMCO  Epoxy 


The  shape  of  the  filament  cross  section  could  be  used 
to  good  advantage  by  a  designer  with  detailed  knowledge  of 


x 

y 

1.345 

.555 

1.21 

.326 

1.15 

.398 
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the  expected  loading  of  the  structure  under  design.   Such 
large  variations  in  some  composite  properties  due  to 
relatively  small  departures  from  a  circular  cross  section 
also  make  it  extremely  dangerous  to  load  a  "circular" 
filament  near  its  expected  transverse  or  shear  strength. 
Small  manufacturing  errors  can  lead  to  large  strength 
discrepancies^ 
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V.   COMPARISON  OF  MATERIAL  COMBINATIONS 

The  preceeding  sections  have  investigated  the  effect 
of  various  constituent  properties  on  the  properties  of 
the  composite.   The  goal  was  to  be  able  to  predict  the 
comparative  properties  of  various  actual  material  combina- 
tions.  In  this  section  several  functional  material  com- 
binations are  analyzed.   Theoretical  mechanical  properties 
are  compared  and  the  results  are  correlated  with  the  guide 
lines  developed  in  the  preceeding  sections. 

Forty,  fifty,  and  sixty  percent  filament  volume  models 
are  compared  for  the  following  combinations  of  materials, 

i)   MODMOR  II/NARMCO  Epoxy 
ii)   MODMOR  11/4  617  Epoxy 
iii)   MODMOR  1/4617  Epoxy 
iv)   Boron/6061  Aluminum 
For  MODMOR  I/NARMCO  Epoxy  and  E-glass/Epoxy ,  3  0%  and  7  0% 
filament  volume  models  were  also  considered.   Single  models 
of  S-glass/NARMCO  Epoxy,  THORNEL  25/4617  Epoxy,  THORNEL 
40/4617  Epoxy,  and  Boron/2024  Aluminum  were  also  investigated 

Primary  measures  of  structural  efficiency  of  a  material 
are  the  specific  strength  and  specific  moduli  of  the  material 
These  are  defined  in  section  I-B,  however  their  importance 
necessitate  further  explanation. 

In  general,  design  of  a  particular  structure  requires 
certain  material  properties  of  strength  and/or  stiffness. 
In  modern  technology  the  weight  of  a  structure  has  gained 
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equal  importance  with  mechanical  properties.   It  is 
significant  to  note  that  all  common  structural  materials, 
when  considered  on  a  stiffness  to  weight  basis,  are  remark- 
ably similar.   Wood,  steel,  aluminum,  and  titanium  fall 
within  the  range  of  92.0  x  10  6  psi/( lb/in3)  to 
105.5  x  106  psi/(lb/in3)  for  specific  Young's  modulus. 
Similarly  for -specif ic  ultimate  strength  these  common 
materials  range  from  4.60  x  105  to  8.33  x  105  psi/ (lb/in3 ) . 
This  data  is  summarized  in  Table  XVI.   Thus  for  a  composite 
to  be  valuable  from  a  structural  standpoint  it  must  exceed 
these  values  for  specific  mechanical  properties. 


TABLE  XVI.   SPECIFIC  ULTIMATE  STRENGTH  AND  SPECIFIC 
YOUNG'S  MODULUS  FOR  COMMON  STRUCTURAL  MATERIALS 


Specific  Ultimate  Strength   Specific  Modulus 
3)xl0"5         psi/(lb/in3)xl0~6 

105.5 

104.2 

97.2 
101.6 

92.0 


Figure  (32)  shows  the  transverse  stiffness,  E  ,  and 
the  specific  transverse  stiffness,  E  /p,  where  p  is  the 
composite  specific  weight.   These  results  are  for  the  50% 
filament  volume  models  considered.   In  the  preceeding 
sections  transverse  stiffness  was  found  to  depend  primarily 
on  the  matrix  modulus  of  elasticity.   The  results  of  Fig. 
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Material    P 

si/  (lfi/i: 

Steel 

5.16 

Aluminum 

6.64 

Titanium 

8.33 

Magnesium 

5.00 

Wood 

4.60 

(32a)  show  that  the  transverse  stiffness  of  materials  with 
similar  matrix  properties  showed  little  variation.   For 
example,  those  composites  with  a  matrix  of  NARMCO  epoxy 
had  comparable  values  of  transverse  stiffness  regardless 
of  filament  properties.   The  aluminum  matrix  was  so  much 
stiffer  than  any  other  matrix  material  studied  that  the 
performance  of  Boron/Aluminum  far  exceeds  that  of  the 
other  composites. 

Figure  (32b)  shows  the  specific  transverse  stiffness. 
Note  that  the  smallest  value  of  composite  specific  transverse 
stiffness  is  approximately  80%  greater  than  that  normally 
associated  with  structural  materials. 

Figure  (33)  shows  the  composite  transverse  strength,  S  , 
and  specific  transverse  strength,  S_/p.   The  constituent 
property  that  showed  the  most  influence  on  transverse 
strength  was  the  matrix  modulus  of  elasticity.   Figure  (33a) 
reinforces  this  judgement  by  showing  negligible  variation 
in  transverse  strength  among  those  composites  utilizing 
identical  matrix  materials.   Progressively  stiffer  matrix 
materials  produced  greater  transverse  yield  strength  regard- 
less of  the  filament  material. 

Figure  (33b)  shows  the  specific  transverse  yield 
strength  for  the  composites  considered.   Note  that  not 
all  of  the  material  combinations  exceed  the  stated  value 
of  8.33  x  105in.  for  common  structural  materials.   One 
reason  for  this  is  that  the  quoted  values  for  common 
structural  materials  are  ultimate  strengths  rather  than 
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Figure  (32) .   Transverse  Stiffness  Comparison. 
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Figure  (33) .   Transverse  Strength  Comparison 
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yield  strengths.   Another  reason  is  that  the  composite 
lamina  are  designed  such  that  their  primary  strength  is 
in  the  longitudinal  direction  whereas  most  standard  struc- 
tural materials  are  considered  isotropic. 

Figure  (34)  is  a  plot  of  longitudinal  stiffness  and 
specific  longitudinal  stiffness.  It  has  been  shown  in 
this  study"  that  composite  longitudinal  stiffness  can  be 
accurately  predicted  by  a  law  of  mixtures.  For  the  50% 
case  of  Fig.  (34a)  the  composites  should  be  in  order  of 
increasing  sum  of  matrix  and  filament  moduli  of  elasticity. 

This  is  seen  to  be  true. 

Figure  (34b)  shows  the  specific  longitudinal  stiffness. 
Note  that  the  graphite/epoxy  combinations  are  considerably 
more  efficient  than  the  Boron/Aluminum  combinations  despite 
the  large  stiffness  advantage  of  aluminum  over  epoxy. 

Figure  (35)  shows  the  longitudinal  yield  strength  and 
specific  longitudinal  yield  strength  of  the  various 

composites . 

The  relative  longitudinal  yield  strength  of  a  composite 
is  somewhat  more  difficult  to  predict  than  are  the  other 
properties.   Increasing  the  filament  modulus  or  the  matrix 
yield  strength  has  been  shown  to  increase  the  composite 
longitudinal  yield  strength.   Increasing  the  matrix  modulus 
of  elasticity  decreases  longitudinal  yield  strength.   There 
is  no  set  method  however  for  predicting  the  combined  effect 
of  such  changes. 
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Figure  (34) .   Longitudinal  Stiffness  Comparison, 
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Figure  (3  5) .   Longitudinal  Strength  Comparison, 
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The  MODMOR  II  (graphite) /Epoxy  models  indicate  a 
case  where  the  matrix  yield  strength  was  of  greater 
significance  than  the  matrix  modulus  increase.   The  4617 
Epoxy  had  a  greater  modulus,  which  should  decrease  the 
yield  strength,  and  a  greater  matrix  yield  strength  which 
should  increase  the  composite  longitudinal  strength.   The 
net  result  was- a  marked  increase  in  composite  yield 
strength.   Considering  the  variations  individually,  as 
in  Fig.  (36) ,  this  indicates  that  the  increase  in  matrix 

yield  strength,  say  from  YQ  *  to  Y   2,  increased  the  composite 

i     m  2    m 
yield  strength  from  S    to   ST .   This  is  shown  in  Fig.  (36) . 

The  decrease  in  matrix  modulus  results  in  a  composite  stress- 
strain  curve  with  a  smaller  slope.   For  the  same  matrix 

yield  strength,  Y  2,    the  composite  yield  is  S3.   The  results 

m 
of  the  analysis  shows  that 


S1  <  S3  <  S2 


(V-l) 


Combined  Effect  of  Matrix  Modulus  and 
Matrix  Yield  (MODMOR  II/Epoxy) 
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The  two  boron  analyses  represent  the  opposite  case. 
The  change  in  matrix  modulus  prevails  over  the  change  in 
matrix  yield  strength  as  the  epoxy  composite  has  a  yield 
strength  considerably  larger  than  that  of  the  Boron/Aluminum 
composite. 

It  should  be  noted  that  the  longitudinal  yield  strength 
of  the  graphite  composite  MODMOR  1/4617  Epoxy  is  considerably 
less  than  that  of  the  MODMOR  11/4  617  Epoxy.   This  is  signif- 
icant in  that  the  Young's  modulus  of  MODMOR  I  is  50%  greater 
than  that  of  MODMOR  II.   This  apparently  contradicts  the 
guidelines  presented  in  this  study.   The  reason  for  this  is 
that  the  conclusions  of  this  investigation  are  based  on  the 
assumption  that  initial  yield  occurs  in  the  matrix  rather 
than  the  filament.   Due  to  the  relatively  low  yield  strength 
of  MODMOR  I  initial  yield  actually  occurs  in  the  filament. 
This  is  predicted  by  considering  the  stress  concentration 
factors  as  outlined  in  section  III-B-6.   Explicitly,  for 
MODMOR  1/4  617  Epoxy, 

Filament  yield  strength   2.5  x  105  psi 

Matrix  yield  strength      .08  xlO5  psi 

Filament  SCFZ  2.106 

Matrix  SCFZ  .0524 


t   =  2.106   =  40  2  (V-2) 

SCFZm     .0524    qu'z  ^V  z; 


Yo 


f   2  5 

=  31.25  (V-3) 


YQ     .08 

m 
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SCFZ^  Yrt 

f  °f 

> (V-4) 

SCFZ  Y_                         v    ' 

m     °m 
m 

Which  meets  the  criteria  of  section  III-B-6  for  initial 

yield  occurring  in  the  filament. 

Figure  (35b)  represents  the  specific  longitudinal  yield 

strength.   The_ unshaded  area,  of  the  bar  chart  represents 

the  specific  yield  strength  from  the  yield  calculations  of 

this  analysis.   The  shaded  region  of  the  figure  represents 

an  estimate  of  specific  ultimate  strength.   The  ultimate 

strength  was  estimated  by  use  of  the  equation 

S  ,A  =  v.Y«   +  vm    a  (V-5) 

Ult     f  of     m  mult 

As  discussed  in  section  III-B-4  am    is  the  stress  in  the 

ult 

matrix  at  the  ultimate  strain  of  the  filament.   Again  it 
should  be  noted  that  this  relation  assumes  that  the  matrix 
remains  elastic  up  to  the  failure  of  the  filament.   Generally 
this  is  not  true.   Using  the  method  of  Lin,  Salinas,  and 
Ito  [9]    the  theoretical  ultimate  tensile  strength  of  a  50% 
Boron/Aluminum  composite  was  found  to  be  243,000  psi.   This 
compares  with  an  estimate  using  Eq.  (V-5)  of  305,000  psi,  a 
difference  of  greater  than  20%.   While  the  discrepancy  is 
rather  large  the  results  are  considered  to  be  worthy  of 
inclusion  as  a  relative  measure  of  the  range  between  yield 
and  failure. 

For  MODMOR  I,  MODMOR  11/4617  Epoxy  and  THORNEL  25  (all 
graphite  filaments)  yield  actually  occurred  in  the  filament. 
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This  is  the  reason  that  the  specific  ultimate  and  specific 
yield  strength  are  virtually  coincident. 

Also  of  importance  in  a  manufacturing  process  is  the 
cost  of  the  materials.   Table  XVIII  shows  strength  per 
unit  cost  data  for  selected  materials.   Cost  data  was 
obtained  from  Refs.   J2  5J  and  [26j  .   The  extremely  high 
strength  to  cost  ratio  for  glass  filament  serves  to  explain 
why  glass  composites  are  in  such  wide  use  despite  their 
relatively  low  transverse  yield  strength.   The  present 
interest  in  graphite  filaments  can  certainly  be  justified 
by  the  strength/cost  increase  anticipated  in  the  near  future 

TABLE  XVII.     COMPARATIVE  STRENGTH  PER  UNIT  COST 
Epoxy  Matrix       v^  =  0.5 
Filament  Strength/Unit  Cost  (in-lb/dollarxlO-1* ) 

Boron  .98 

Graphite  2.01 

Graphite  (1974)                 4.02 

Glass  42.48 

Numerous  other  factors  also  enter  into  a  material 
selection  process.   Ease  of  handling  and  fabrication, 
special  processing  requirements,  and  specific  desirable 
properties  such  as  thermal  expansion  may  be  such  factors. 
Consideration  of  these  factors  is  beyond  the  scope  of  this 
investigation . 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

In  this  investigation  a  finite  element  analysis  of  the 
microscopic  and  macroscopic  stresses  of  a  composite  material 
has  been  carried  out.   The  stresses  were  used  to  calculate 
the  initial  yield  surface  for  unidirectional  composite 
materials.   The  results  of  these  analyses  were  then  used 
to  determine  the  mechanical  properties  of  a  laminate  of 
composite  material.   Parametric  studies  of  the  effect  of 
various  constituent  properties  on  the  composite  mechanical 
properties  were  then  carried  out.   Table  XVIII  summarizes 
the  results  of  these  parameter  variations. 

These  results  are  based  on  the  assumption  that  initial 
yield  takes  place  in  the  matrix  material.   The  effects  noted 
in  this  table  are  those  associated  with  a  change  of  the 
specified  constituent  property  only,  that  is,  with  all 
other  constituent  properties  held  constant.   Caution  must 
be  exercised  in  attempting  to  predict  the  effects  of  changes 
in  more  than  one  constituent  property. 

Other  results  included: 

i)   Composite  Transverse  Poisson's  Ratio,  v  .   This 
value  exceeded  that  of  either  constituent.   By  adjusting 
the  ratio  of  filament  Young's  modulus  to  matrix  Young's 
modulus,  E  /Ef,  the  designer  may  achieve  extreme  values  of 
transverse  Poisson's  ratio. 
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ii)   Longitudinal  Initial  Yield  Strength,  SL.   For 
initial  yield  occurring  in  the  matrix  material  the  longi- 
tudinal yield  strength  may  be. predicted  by 

S   =  V  Of  +  vmYQ 
Jj  m 

where  Y    is  the  yield  strength  of  the  matrix  material, 

m 
a£  is  the  stress  in  the  filament  at  the  yield  strain  of 

the  matrix,  and  v^  and  v   are  filament  and  matrix  volume 
fractions,  respectively. 

iii)   Initial  Yield  Location.   The  computer  program 
utilized  in  this  investigation  gives  results  which  include 
the  location  of  the  yield  point  for  any  macroscopic  stress 
state.   For  load  case  III  it  may  be  determined  whether  the 
composite  fails  in  the  matrix  or  the  filament  by  consider- 
ing the  stress  concentration  factor  in  the  longitudinal 
direction,  SCFZ.   If  the  ratio  of  SCFZ  in  the  filament  to 
that  in  the  matrix  is  less  than  the  ratio  of  filament  yield 
strength  to  matrix  yield  strength  then  yield  will  occur  in 
the  matrix.   That  is,  if 

SCFZf/SCFZ   <  YQ  /YQ 

f   m 

yield  will  occur  in  the  matrix.   This  determines  the  appli- 
cability of  the  longitudinal  yield  equation  noted  above  and 
does  not  require  the  application  of  any  particular  yield 
criteria. 

iv)  Filament  Cross  Section  Effects.  Variations  from  a 
circular  filament  cross  section  have  significant  influence 
on  composite  properties.   Shear  modulus  and  transverse 
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properties  change  considerably.   Due  to  the  distribution 
of  microscopic  stresses  on  the  block  boundaries  an  el- 
liptical cross  section  has  greater  transverse  strength  in 
the  direction  of  the  minor  axis. 

B .   RECOMMENDATIONS 

i)  The  logical  extension  of  this  investigation  is 
consideration  of  inelastic  loading.  The  analysis  used 
in  this  study  is  a  special  case  of  the  elastic-plastic 
analysis  of  Lin,  Salinas,  and  I  to.  [9]  [Lu] 

ii)   Further  effort  seems  warranted  in  the  evaluation 
of  filament  cross  section  effects.   Through  use  of  cross 
sectional  shapes  and  filament  distribution  it  should  be 
possible  to  exert  control  over  transverse  and  shear 
properties. 

iii)   The  wide  variation  in  the  literature  of  theoretical 
transverse  Poisson's  ratio  certainly  suggests  an  urgent 
need  for  accurate  experimental  results  as  well  as  a 
comprehensive  survey  of  analytical  techniques. 
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APPENDIX  A:   YIELD  LOCUS  CALCULATIONS 

The  analytical  method  employed  in  this  investigation  8 
requires  that  the  macroscopic  stress  state  of  the  composite 
be  one  of  plane  stress.   It  is  also  required  that  the  rec- 
tangular form  of  the  basic  block,  Figs.  (2)  and  (3),  be 
maintained  under  deformation.   To  ensure  these  results 
four  distinct  loading  cases  are  superposed.   These  load 
cases  and  their  boundary  conditions  are  shown  in  Fig.  (4). 

Load  cases  I  and  II  are  plane  strain  states  with  e      =0. 

r  z 

Load  case  III  is  a  generalized  plane  strain  case  with 
E   =  1.0.   Load  cases  I,  II,  and  III,  are  superposed  to 
give  Sy  =  0.   Load  case  IV  is  a  state  of  longitudinal  shear. 

All  four  cases  have  S    =  S    =0.   Rectangular  shapes  of 

xy     yz  3         c 

the  basic  block  result  from  all  load  cases. 

Using  finite  element  methods  with  linear  strain  triangles 
[l2J  the  four  loading  cases  are  individually  analyzed  to 
obtain  the  microstresses  and  their  associated  macrostresses. 
The  microstresses  are  designated  a-   •  with  a  subscript  denoting 
the  stress  component  and  a  superscript  for  the  load  case. 
For  example,  for  load  case  I  the  microstresses  obtained  are 

ax'  axy'  ay'  az  (A_1) 

Similar  results  are  obtained  for  problems  II  and  III.   The 

macrostresses  for  load  case  I  are  designated 

Sj,  S?,  S1  (A-2) 

x '   z '      y 

and  similarly  for  load  cases  II  and  III. 
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Load  case  IV  yields  microstress  results  for 


a1 


v 


xz 


and 


<jIV 

yz 


(A-3) 


and  macrostresses 


SIV  and  SIV 
xz      yz 


(A-4) 


In  the  notation  of  Lin,  -Salinas  and  Ito  [_8J  the  micro- 
stresses  of  a  combined  loading  state  may  be  related  to  the 
macrostresses  as 


and 


{a} 
4x1 


{a} 
2x1 


[a]  {s} 

4x3  3x1 


[a]  {S} 
2x2  2x1 


(A-5) 


(A-6) 


Where 


xa '  x  y  xy  z 


{S  }   =  <  S  S  S   > 
x      s  x  y  z 


(A-7) 
(A- 8) 


{a}   -<  oxzoyz> 
-   T 

<s>   =  < sxz  syz 


(A-9) 
(A-10) 


[Aj  - 


All  A12  A13 
A21  A22  A23 
A31  A32  A33 
A41  A4  2  A4  3 


(A-ll) 
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M  - 


All  A12 


(A-12) 


A21  A22 

The  matrix  A  may  be  related  to  the  results  of  load 
cases  I,  II  and  III  as 


I   II   III 

ax  ax   ax 


I  II  III 

ay  ay  aY 


I   ii  in 

axy  axy  axy 


I   II   III 

°z  az   az 


An  A12  A13 


A21  A22  A23 


A31  A32  A33 


A41  A42  A43 


Sx  *lZ    <1T 


s1  s11  s111 
sy  sy   y 


(A-13) 


For  convenience  this  may  be  written 


[T]  =  [A]  W 
4x3   4x3  3x3 

The  matrix  A  can  then  be  evaluated  as 


w  =  mm:1 

4x3    4x3  3x3 


(A-14) 


(A-15) 


Thus  the  combined  microstresses ,  [a]  ,  due  to  any  uniform 
transverse  or  longitudinal  deformation  may  be  evaluated  by 
substituting  the  results  of  Eq.  (A-15)  into  Eq .  (A-5)  as 


{o}      [T]   [C1  <S> 
4x1      4x3   3x3   3x1 


(A-16) 


in  a  similar  manner  for  load  state  IV  and  an  analogous 
load  case  V  associated  with  S^  shear  loading,  the  matrix 

A  may  be  evaluated  as 
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[a]  - 

IV 

°xz 

V 
axz 

2x2 

IV 

V 
ayz 

.IV 
*xz 


;V 
yz 


(A-17) 


This  may   be  written   as 

pa  =  m  w 

2x2   2x2   2x2 


(A-18) 


From  Eqs.  (A-6)  and  (A-18)  the  microstresses ,  {a},  may  be 

related  to  any  specified  shear  loading  as 

{a}  =[t]   W"1  <S}  (A-19) 

2x1   2x2   2x2    2x1 

The  derivation  of  Eqs.  (A-16)  and  (A-19)  is  independent 
of  the  load  state.  For  the  particular  case  of  plane  macro- 
scopic stress  treated  in  this  investigation  the  macroscopic 
stress  vectors  become 

{S}T  =  <sx  0  sz> 


{s>T  =  <sxz  0  > 


(A-20) 

Equations  (A-16)  and  (A-19)  may  be  evaluated  for  the  point 
by  point  microstresses  associated  with  a  plane  macroscopic 
load  condition  as  specified  by  Eq.  (A-20) . 

Von  Mises  yield  criterion  is  assumed  to  govern  initial 
yield.   This  requires  that   1 

-  k2  +  a 

1      i-    X        ' 


J  =  a   2  +  a_2  +  a  __2 


xz 


yz 


xy 


"  °x°y  "  V«  "  a^x]  =  k     Y02 


(A-21) 


Y   is  the  elastic  limit  of  the  material  for  the  tensile 
o 

stress  state. 
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Equations    (A-16) ,     (A-19) ,     (A-20) ,    and    (A-21)    may   be 
combined   to   form  a   yield   equation    for   plain  macroscopic 
loading   in   terms   of   the  macrostresses   and   the  point   depen- 
dent  components   of    the      A     matrix.      This    equation   is 

[A33+    3(A13+A23+A43_A13A23"A23A43"A43A13)J    Sz 
-    +     [2A31A33+   |    ^AHA13+A21A23+A41A43) 

"    3    (AHA23+A13A21+A21A43+A23A41+A41A13+A43A11)]     SzSx 


[A31+   I 


(a£1+a|1+a|1-a11a21-a21a41-a41a11 


s.2. 


+  (A^  A|x)  S2Z  -  1  Y2  =  0  (A-22) 

For  given  values  of  S   and  Sv_  Eq.  (A-22)  may  be  solved 

X        X  z 

for  the  value  of  S   necessary  to  cause  yield.   Each  point 

within  the  fundamental  block  will  give  a  different  value  of 

S   dependent  on  its  microstress  state.   Each  combination  of 

Sv  and  S„„  yields  two  values  of  S_  at  each  point,  a  maximum 

and  a  minimum.   The  minimum  of  the  maximum  values  and  the 

maximum  of  the  minimum  values  determine  two  yield  points  for 

the  given  values  of  Sv  and  S   .   Continuing  in  this  manner 
3  x      xz  3 

the  initial  yield  surface  for  a  composite  subjected  to  a 
macroscopic  plane  stress  state  is  calculated. 
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APPENDIX  B:   CONSTITUTIVE  EQUATIONS 
FOR  A  HOOKEAN  FIBROUS  COMPOSITE  MATERIAL 


The  well  known  equations  of  elasticity  relating  stress 
components  and  strain  components  in  terms  of  two  independent 
elastic  constants  represents  a  special  case  of  generalized 
Hooke ' s  Law  for  isotropic  materials.   Here  we  derive  the 
Hookean  stress-strain  relations  for  a  lamina  of  an  anisotropic 
fiber  reinforced  composite  material. 

The  generalized  Hookean  stress-strain  law  relating  stress 
and  strain  may  be  written  in  matrix  form  as 


a    a    a    a    a    a 
11   12   13   14   15   16 

a21  a22  a23  a24  a25  a26 


X 

ey 

E2 

Gxy 

"■" 

eyZ 

xz 

a„,  a_^  a„_  a„  .  a„_  a^^ 
31   32   33   34   35   36 

a    a    a    a    a    a 
41   42   43   44   45   46 

a   a   a   a   a   a 
51   52   53   54   55   56 

a    a    a    a    a    a 
61   62   63   64   65   66 


"s 

X 

sy 

s 

z 

S 

xy 

S 

yz 

S 

xz 

(B-l) 


In  a  more  compact  notation  this  may  be  written 

{e}  =  [a]  {S} 
6x1   6x6   6x1 


(B-2) 


The  £• .  terms  are  engineering  strain  components  and  the  S. . 
are  stress  components. 

Thus  the  generalized  stress-strain  relations  contain  36 
unknown  Hookean  constants, 

a±j  (i,j  =  1,2,3,4,5,6)  (B-3) 

Consideration  of  the  strain-energy  density  function  J27J 
establishes  symmetry  of  the  matrix  of  Hookean  constants 


aij  "  aji 


(B-4) 
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and  reduces  the  maximum  number  of  independent  elastic 
constants  to  21.   Further  reductions  in  the  number  of 
elastic  constants  depends  on  the  number  of  planes  of 
symmetry  present  in  the  material. 

When  the  elastic  constants  at  a  point  have  the  same 
values  for  every  pair  of  co-ordinate  systems  which  are 
mirror  images  of  each  other  in  a  certain  plane  that  plane 
is  called  a  plane  of  elastic  symmetry.   This  means  that  if 
a  certain  plane,  say  the  (x-,  ^2)  plane  in  co-ordinate 
system  (xpXjfXj  is  a  plane  of  elastic  symmetry  then  the 
reflection  of  this  system,  (x',x',x'),  where 

x'  =  x 


1  1 

2  =  x2 


x3  =  -x3  (B-5) 

must  have  the  same  elastic  constants.   If  you  could  actually 
"see"  the  elastic  properties  this  means  that  the  properties 
would  appear  identical  at  every  point  regardless  of  which 
side  of  the  plane  you  were  viewing.   Alternatively  Sines  [2$] 
suggests  that  if  you  could  actually  see  these  properties 
they  would  look  the  same  even  if  the  observer  stood  on  his 
head,  that  is  rotated  180   about  his  line  of  sight.   This 
invariance  is  utilized  to  evaluate  the  elastic  constants. 

Consider  the  reference  system  (x-,  ,x„,xJ  of  Fig.  (B-l)  . 
If  the  plane  (x^,x?)  is  a  plane  of  elastic  symmetry  a 
reflection  through  the  x^  axis  may  be  carried  out  without 
effecting  the  elastic  constants.   This  new  reference  system 
is  designated  by  the  primed  axes  of  Fig.  (B-l). 
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x2,x3 


/ 


/ 


/ 


X.  .X 


]/*l 


Figure  (B-l) .   Reflection  of  a  Plane  of  Symmetry  (x1/x2). 
The  general  equations  governing  the  transformation  of 


axes  are 

-  ,      0  2„      ,_   2„   L„   2 


ll   ex+ml  cy+ni  ez+£lmlexy+minleyz+ni£iezx 
ey  "  £22ex+m22£y+n22£z+£2m2£xy+m2n2£yz+n2£2ezx 


.32£x+m32ey+n32 


In    ev+  n0ze„+n^^£z+£3m3ex   +m3n3e      +n3il3£zx 


£ 


£     s  £ 


Gxy  =    2Wx   +2m1m2£y   +  2nXn2  %+ <  V2+IV2>  Sty 


+    (m1n2+n1m2)e      +(n1£2+£1n2)e 


zx 


r-yZ  =  2£2£3£  +2m2m3e  +2n2n3ez+ (£2m3+m2£3)  exy 
+  (m2n3+n2m3)eyz+(n2£3+£2n3)£zx 

Ezx  ■  2^3^iex+2m3mi£y+2n3nl£z+(£3ml+illm3)exy 

+  (m3n1+n3m1) eyz+ (n3£1+£3n1)  ezx 
for  strain  components  and 


(B-6) 
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Sx  =    4Sx+mISy+niSz+2£lmlsxy+2minlSyz+2nl£lSzx 


Sy   -    £2sx+ni|sy+n|sz+2£2m2Sxy+2m2n2Syz+2n2^2Szx 

Sz   =    *3V-m3^y+n32V2£3m3Sxy+2m3n3Syz+2n3£3szx 

Sxy   =    £l£2Sx+mlm2Sy+nln2Sz+(V2+ml£2)Sxy 

+    (m1n2+nim2)Syz+(n1£2+£1n2)Szx  (B-7) 

Syz    "    £2£3Sx+m2m3Sy+m2n3Sz+(^2m34-m2^3)Sxy 

+    (m2n3+n2m3)Syz+(n2£3+£2n3)Szx 

szx   =    ^3^1sx+m3mlSy+n3nlsz+(£3ml+m3£l)sxy 

+    (m3n1+n3m1)Syz+(n3S,1+£3n1)Szx 

for   stress   components.      The   terms 

lif   m±,    ni         (i=  1,2,3)  (B-8) 

are   direction   cosines    defined   by    the    following    table: 

12         3 
1' 

2' 

3' 

This  indicates  that  1^  cos  (x|  ,x-|_)  etc.   For  the  reflection 
of  Fig.  (B-l)  the  direction  cosines  are: 

12    3 

1" 

2' 

3' 


h 

mi 

nl 

H 

m2 

n2 

i3 

m3 

n3 

1 

0 

0 

0 

1 

0 

0 

0 

-1 
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Using  these  direction  cosines  to  evaluate  the  transformation 
relations  Eqs.  (B-6)  and  (B-7)  yields 

sx  =  sx         ex  =  ex 

Sy  =  sy       ey  =  ey 

SZ  =  SZ  £Z  =  £2 

Sxy  =~Sxy       £xy  =  ~£xy  (B-9) 

<5  •    =—  Q  £     as  —  £ 

yz   yz      yz    yz 

Szx  =  Szx       £zx  =  ezx 

The  stress-strain  relation  in  the  new  co-ordinate  system 

may  be  written 

T  T 

<ex'  Ey  ez'  ex'y  ey'z  ez'x  >  =  [a]  <Sx  Sy  Sz  S   -Syz  -Szx>    (B-10) 

6x6 

The  elastic  constants  are  then  evaluated  by  equating  the 

transformed  strains  to  the  original  strains  in  accordance 

with  Eq.  (B-9).   Consider  for  example  ej  =£„•   In  detail 

this  requires  that 

allsx+a12Sy+a13Sz"a14Sxy"a15syz+a16Sxz  = 

a11Sx+a12Sy+a13Sz+al4Sxy+a15Syz+a16Szx       (B-ll) 

The  constants  are  invariant  with  transformation,  requiring 

that 

"a15  Syz  =  a15  syz 

-a14  Sxy  =  a14  Sxy  ^B"12) 

For  a  non  zero  stress  condition  this  can  be  true  only  if 

a15  =  a14  =  °  (B-13) 

Carrying  this  out  for  the  other  relations  of  Eq.  (B-9) 
shows  that 

a15=a14=a24=a25=a34=a3  5=a46=a56=0  (B-14) 
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Thus  a  single  plane  of  elastic  symmetry  reduces  the 
number  of  independent  elastic  constants  to  13. 

If  a  material  has  three  orthogonal  planes  of  symmetry 
a  reflection  of  a  second  plane  may  be  carried  out  without 
changing  the  Hookean  coefficients.   Reflecting  the  (x2,x-J 
plane  through  the  x,  axis  results  in  the  transformed  axes 
of  Fig.  (B-2) .   In  this  case  the 


x3,x3' 


Figure  (B-2).   Reflection  of  a  Plane  of  Symmetry  (x2,x  ) 


direction  cosines  are  given  by 


2    3 


1" 

2' 
3' 


-1 

0 

0 

0 

1 

0 

0 

0 

1 

The  transformation  relations  Eqs.  (B-6)  and  (B-7)  give 


125 


s 

X 

— 

S 

X 

s; 

= 

S 

y 

Sz 

= 

Sz 

■ 

-s 

xy 

X 

i 
yz 

= 

syZ 

■ 

_ 

-s 

zx 

z 

x    ; 
i 

£  =  £ 

y   : 

£  =  £ 


Z 

-£ 


xy 
i 

yz 
zx 


xy 


=  £ 


yz 

-£ 


ZX 


(B-15) 


Using  these  relations  it  may  be  shown  that 

a26  =  a16  =  a36  =  a45  =  ° 
The  matrix  of  elastic  constants  may  now  be  written 


(B-16) 


M. 


6x6 


all  a12 

a13 

0 

0 

0 

a22 

a23 

0 

0 

0 

a33 

0 

0 

0 

a44 

0 

0 

symmetric 

a55 

0 

(B-17) 


a66 

where  the  subscript  o  indicates  that  this  matrix  is  defined 
for  an  orthotropic  material.   There  are  now  nine  independent 
elastic  constants. 

Further  reductions  in  the  number  of  independent  elastic 
constants  depend  upon  the  number  of  planes  of  isotropy,  if 
any,  present  in  the  material.   A  plane  of  isotropy  is  one 
in  which  the  elastic  constants  are  invariant  with  respect 
to  direction  within  the  plane. 

Assume  that  a  plane  of  isotropy  does  exist,  the  (x-,  ,^2^ 
plane  of  Fig.  (B-3) .   This  indicates  that  the  Hookean 
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coefficients  are  independent  of  the  orientation  of  the 
(x  ,x9)  axes  within  this  plane.   An  infinite  number  of 
orientations  could  be  chosen  within  this  plane  to  assist  in 
evaluating  the  remaining  unknowns .   A  45  degree  rotation, 
shown  in  Fig.  (B-3) ,  was  selected  for  ease  of  calculation. 


x„  .x 


3'  3 


Figure  (B-3).   45  Degree  Rotation  of  Plane  (x^,x2) 

The  direction  cosines  for  such  a  case  are: 

12    3 

1' 


2\ 


3' 

Again  using  the  transformation  relations  Eqs.  (B-6)  and 
(B-7)  one  obtains 

ex  =  h   ex  +  *  ey  +  %   exy 
ey  "  **    ex  +  h    £y  "  h    exy 


1 

r2 

1 

/2 

0 

i 

/2 

1 

7? 

0 

0 

L         . 

0 

1 

(B-18) 
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^xy         ex         ey 


yz        /2      yz     /2    zx 

zx      /j    fcyz    T/2      zx 


S '    =   -  S      +   i  S      +S 
x         2      x        7      y  xy 

Sy   =  I  Sx   +  I  Sy   "   Sxy 


s'   =  s 

z  z 

S'       =   I  S      -is 
xy        2      x        2      y 

;' 
yz 


S'      =  -is       +is 

/J     yz     /2     zx 


(B-18) 


S'      =  i  s       +  1  S 
zx      /J     yz      /?      zx 

Using  these  relations  it  may  be  shown  that 


a55  "  a66 


all  "  a22 


a13  "  a23 

a44  -  2(a11  -  a12) 


(B-19) 


A  material  exhibiting  one  such  plane  of  isotropy  is 
called  a  transversely  isotropic  material.   There  are  five 
remaining  independent  elastic  constants  and  the  matrix  of 
elastic  constants  may  be  written 


[• 


TI 


6x6 


all  al2  al3 

0 

0 

0 

all  ai3 

0 

0 

0 

= 

a33 

0 

0 

0 

2(a11-a 

12) 

0 

0 

symmetric 

a« 

0 
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55 


(B-20) 


where  the  subscript  TI  signifies  transversely  isotropic. 
The  stress-strain  relations  for  this  case  may  be  written 

U}    '    [ATlJ   {S}  (B-21) 

6x1     6x6   6x1 

One  of  the  assumptions  of  this  study  is  that  a  composite 
material  may  be  considered  to  be  homogeneous  in  a  macro- 
scopic sense.   This  assumption  allows  the  composite  to  be 
considered  as  a  transversely  isotropic  material,  and  Eqs. 
(B-21)  are  the  constitutive  equations  for  a  composite 
material.   Introducing  the  engineering  constants  these 
equations  become 


e„  =  — L  (Sv  -  vms  )  -  J±      S 
E 


x    „   v  x     t  y     e~    z 


T 

S, 
L 


y   et  v  y   t  x'   E 


Z*  =  '   e7  (sx  +  sy)  +  e7  Sz  ^b"22) 

L  Xj 

£    =     c 

xy   G^  xy 

yz     gt    yz 

£    =   l    S 
J-l 

ET 
GT     2(1  +  vT) 


E   and  G   are  the  Young's  modulus  and  shear  modulus  in 
the  transverse  direction,  in  the  plane  of  isotropy.   E   and 

Xj 

G   are  corresponding  values  in  the  longitudinal  direction, 

Xj 

perpendicular  to  the  plane  of  isotropy.   Poisson's  ratio 


129 


vT  signifies  the  contraction  in  the  plane  of  isotropy  due 

to  tension  in  the  same  plane.   Poisson's  ratio  v   represents 

jj 

the  contraction  in  the  plane  of  isotropy  due  to  tension  in 
the  longitudinal  direction.   The  transverse  shear  modulus 
is  a  dependent  quantity. 
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APPENDIX  C:   YIELD  LOCI 

The  figures  included  in  this  section  represent  the 

yield  loci  that  were  generated  during  the  course  of  this 

investigation.   Each  plot  is  labelled  with  the  constituent 

materials,  filament  volume  fraction.  v_.  and  S,  and  S„ 
"  f       x      z 

scales.   The  outermost  locus  for  all  cases  corresponds 

to  Sv-   =    0.   The  S,,,,  increment  between  successive  loci  is 
««  xz 

also  noted  on  each  plot. 

Those  plots  that  were  used  in  the  parametric  analysis 

are  labelled  in  terms  of  the  parameter  being  investigated. 

All  other  properties  are  held  constant. 

Note  that  in  some  cases  the  S„  and  S   scales  are  dif- 

x      z 

ferent  from  one  another.   This  occurs  when  the  matrix 
materials  are  quite  soft  in  comparison  to  the  filament 
material.   The  transverse  strength  is  so  much  less  than 
the  longitudinal  strength  that  the  surface  plots  nearly 
as  a  straight  line  if  the  same  scales  are  used  for  both 
axes. 
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APPENDIX  D:   PROGRAM  OVERLAY  STRUCTURE 

The  overlay  root-segment  structure  is  shown  in  Fig. 
(D-l)  where  each  rectangle  represents  a  subroutine.   The 
subroutine  functions  are: 

MAIN  -Remains  in  core  and  controls  calling  sequence 
during  execution. 

BANDEC  -  Solves  banded  matrix  problems.   Used  in  MAIN1 
and  MAIN2. 

MULI  -Multiplies  a  symmetric  banded  matrix  with  a 
vector.   Used  in  MAIN1  and  MAIN2. 

MAINl  -  Calculates  microstresses  and  macrostresses  for 
elastic  loading  of  load  cases  I,  II ,  and  III. 

DAVES  -Computes  element  stiffness  and  assembles  system 
stiffness . 

CALTAU  -Calculates  element  stresses  and  strains. 

REVGNK  -  Revises  system  stiffness  to  accommodate 
boundary  conditions  and  applied  loads. 

MAIN2  -  Calculates  microstresses  and  macrostresses  for 
longitudinal  shear  elastic  loading,  load  case  IV. 

MAIN3  -  Calculates  Initial  Yield  Locus. 

MATINV  -  Matrix  inversion  with  accompanying  solution 
of  linear  equations. 
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APPENDIX  E:   COMPUTER  PROGRAM 
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